{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/var/folders/np/m42416h166j90krvxmfl40940000gp/T/ipykernel_54471/3423083904.py:3: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here.\n",
      "Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations\n",
      "  solTchem = (data[1:,:]).astype(np.float)\n"
     ]
    }
   ],
   "source": [
    "data = np.genfromtxt(\"CMAQ_OH_HNO3.dat\", dtype=str)\n",
    "Header = (data[0,:]).tolist()\n",
    "solTchem = (data[1:,:]).astype(np.float)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "campsol = np.loadtxt(\"CMAQ_OH_HNO3_results.txt\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "A_indx = Header.index('A')\n",
    "B_indx = Header.index('B')\n",
    "C_indx = Header.index('C')\n",
    "t_indx = Header.index('t')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of Cells: 1\n",
      "Number of time iterations: 127.0\n"
     ]
    }
   ],
   "source": [
    "niterT, Nvars = np.shape(solTchem)\n",
    "Nsamples = len(np.where(solTchem[:,0]==-1)[0])\n",
    "print('Number of Cells:',Nsamples)\n",
    "print('Number of time iterations:',niterT/Nsamples)\n",
    "solTchem = solTchem.reshape( int(niterT/Nsamples), Nsamples,Nvars)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7fa87051f7f0>"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABDsElEQVR4nO3dd3hUZfbA8e+ZSYdACKETSOjSpYoUUUCKBRuKDUUF69pd8Wfv7uoqYkNEARvoYsOGIE0EVrr03hJqSCCEkjbz/v54JxAgJBOSySSZ83me+8zMnTv3nqtkzrxdjDEopZQKXA5/B6CUUsq/NBEopVSA00SglFIBThOBUkoFOE0ESikV4IL8HUBhxcTEmLi4OH+HoZRSZcqSJUv2G2Oq5fVemUsEcXFxLF682N9hKKVUmSIi28/0nlYNKaVUgNNEoJRSAU4TgVJKBThNBEopFeA0ESilVIDzWSIQkU9EZJ+IrDrD+yIio0Rkk4isEJF2vopFKaXUmfmyRDAe6JfP+/2Bxp5tOPCBD2NRSil1Bj5LBMaYP4CUfA4ZCHxqrP8BUSJSy1fxLNtxgP5vzyXb5fbVJZRSqkzyZxtBHSAh1+tEz77TiMhwEVksIouTkpLO6mLj529j7e5DrNyZelafV0qp8sqfiUDy2JfnKjnGmDHGmA7GmA7VquU5QrpA/VvawkZIkLaPK6VUbv78VkwEYnO9rgvs8tXFghw277i1ZkgppU7iz0QwBRji6T10HpBqjNntq4s5nTYRZGsmUEqpk/hs0jkRmQj0BGJEJBF4FggGMMaMBn4BBgCbgKPAUF/FAuAUmwhcbl2jWSmlcvNZIjDGXF/A+wa411fXP1VO1ZAmAqWUOlnAtJw6HTlVQ5oIlFIqt4BLBFoiUEqpkwVeIjCaCJRSKreASQRBDnurLpcmAqWUyi1gEoG2ESilVN4CLhFoG4FSSp0s4BKBDihTSqmTBUwiOD7FhDYWK6XUSXw2oKy0OV4i0MZipVQZY4whOT0ZQagaXrXYzx8wiSDIqW0ESqnSK9OVyc7DO0lISyAhLYHEtMTjjzsP7yTdlc6wVsO4v939xX7tgEkEOXMNFdhrKDsTEhdC2h5o1BvCo3wfnFIqIGS5skg4nMCOQzvYfmg7Ow7tYEfaDnYc2sGeo3twmxNtmOFB4dSNrEu9SvXoWqcrdSrW4dzq5/okrsBJBN70GnJlwedXwba59nVoZej9DHS8owQiVEqVBznVOFtTtx7fth3axvZD29l5eOdJX/aVQipRv1J92lZvS71K9agXWY/YyFjqRtalalhVRPJatqX4BUwiyBlQlm+J4PfnbBLo+yrUPhfm/At+fgQOJkDv56CE/qcopUo/l9vFrsO72Jy6mS2pW9hycMvxL/60rLTjx4UHhVMvsh7Nqzanf3x/4irFUa9SPepH1icqLMp/N5BL4CSC420EZ+g+mpoIC96F9kOhyz12303fwC+PwryREFIRLnisZIJVSpUabuNm1+FdbDq46fi2+eBmtqZuJcOVcfy4auHViK8cz4AGA2hQuQFxleNoULkB1SOq45DS3UEzYBJBgSOL1/1sH7vkmhnb4YRL3oTMozDrJahSH1pf6+NIlVL+kpqRyoYDG9hwYAMbD2xk44GNbDq4iaPZR48fU7NCTRpWbkinmp1oFNWI+MrxNIhqQKWQSn6MvGgCJhEcX4/gTN1H1/4IMU0hpvHJ+0Xg8ndsiWHK/VCtGdRq7eNolVK+5DZuEtMSWZeyjvUH1rM+ZT3rD6xnz5E9x4+JCo2icZXGXNHoChpXaUyjqEY0jGpIZEikHyP3jYBJBPmWCI4kw/b50O3BvD8cFAKDxsGHF8BXN8Gdf2hvIqXKiGx3NltTt7ImeQ1rU9ayNnkt6w+s50jWEQCc4iS+cjztqrejaXRTmlZpSpMqTYgJjymxxlp/C5hEICI4HZJ3r6GN08C4oNmlZz5Bxepw7QQY1x+m/AOu/VQbj5UqZVxul/3ST1nD6v2rWZ28mvUp60l3pQO24bZJlSZc2uBSzok+h2ZVm9EoqhGhzlA/R+5fAZMIwJYK8iwR7FkJQeFQq23+J4jtBL2egenPwOKPtVupUn5kjGHv0b2sSFrBqv2rWLl/JWuS1xyvzw8PCuec6HO4psk1NK/anOZVmxNXKQ6nw+nnyEufgEoEQQ4h25VHr6H96yGmETi8aNnv8g/Y+gf89iTUOx9qNC/+QJVSp0nPTmd18mpWJK3g76S/WZG0gqRjSQAEO4JpFt2MyxteTsuYlrSMaalf+oUQUIngjCWCpA321743HA644gP4oCtMvg2Gz4Lg8OINVCnF3iN7WZa0jL/3/c3yfctZl7KObJMNQGxkLB1rdqR1tda0jmlN0+imhDhD/Bxx2RVQiSDY6Ti9jSDzCKTugHY3e3+iitXhyg/g86ttNdGA14s3UKUCjDGGralbWbJvCUv3LmXZvmXsPLwTgDBnGC1jWnJLi1toU60Nbaq3ITos2s8Rly8BlQjyLBHs32gfY5oU7mSNesN598L/3oOGvaBpv+IJUqkA4DZuNh7YyOK9i1mydwlL9i4hJT0FgOiwaNrXaM8NzW6gXQ3bkyfYEezniMu3gEoEQQ45fWTx/g32sVrTwp+w97O2veCHe+DuBRBZo+hBKlUOGWPYdHATC/csZNGeRSzeu5jUjFQA6lSsQ7c63Whfoz3tqrejfqX6AdNts7QIqESQZ4kgaT2IA6IbFP6EQaFw9VgY0xO+vwtu/Ma7BmelAkBCWgJ/7f6Lv3b/xcI9C4//4q9TsQ4Xxl5Ix5od6VCjA7Ur1vZzpCqgEkFQXuMI9q+HKvH2S/1sVG8GfV+Gnx+G/70P599X9ECVKoNSM1JZuGchC3YtYMGuBSQeTgTsHDxdanehc83OdKrViToV6/g5UnWqgEoEeZYIkjcXvn3gVB1ug80z7eyl8d2hVpuinU+pMsDldrEqeRXzds5j/q75rNy/ErdxExEUQadanbip+U10qdWF+MrxWtVTygVUIghyOE4fR5C2B+qdV7QT58xH9MH5MPl2uHMOhFQo2jmVKoVS0lOYt3Mec3fOZcGuBRzMOIggtIxpyR2t7uD82ufTulprbdwtYwIqEZw2xYQrG46lQIXqRT95RDRcNQYmXA6//hMGvlf0cyrlZ8YYNhzYwOyE2fyx8w9WJq3EYIgOi6ZH3R50rd2V82ufX2rm1VdnJ6ASQbBTyMo9++jRZPtYIaZ4LhDfA7o/AnPfgAYXQqtriue8SpWgLFcWi/YsYmbCTOYkzjk+I2fLqi25u83d9Kjbg3OqnlPq59hX3guoRHBaieDIPvtYoVrxXaTnE3aVsx8ftKucVW1YfOdWykcOZx5m7s65zNwxk7k753Ik6wjhQeF0qdWFe9rcQ/e63YkJL6YfTKrU8WkiEJF+wNuAExhrjHntlPcrA58D9TyxvGGMGeereIKcDrJytxEcsfOUULEYqoZyOINsl9LR3ewUFLdPO/seSUr5UEp6CrMTZjN9+3T+2v0XWe4sosOi6RfXjwtjL6Rzrc6EBYX5O0xVAnyWCETECbwH9AESgUUiMsUYsybXYfcCa4wxl4lINWC9iHxhjMn0RUxBDiEzO3ci2G8fi7NEABBVz7YRfHWT7UnU79XiPb9SZ2n/sf3M2D6DadunsXjvYtzGTZ2Kdbih2Q30qt+L1jGtdaK2AOTLEkEnYJMxZguAiEwCBgK5E4EBIsX2LasIpADZvgooyOngSKbrxI7DPqgaynHOZdDpTju2IK4bNLuk+K+hlBdS0lP4ffvvTN02lSV7l+A2buIqxXF7y9vpU78PzaKbaffOAOfLRFAHSMj1OhHofMox7wJTgF1AJHCdMea0eaJFZDgwHKBevXpnHVDwqdNQH0kCRzCEVT7rc+br4hch4S/4/m64c65d81ipEpCWmcaMHTP4deuv/LX7L1zGRVylOIa1GkbfuL40imqkX/7qOF8mgrz+lZ06B3RfYDlwEdAQmC4ic40xh076kDFjgDEAHTp0OMOiwwULcgrZuXsNHdlvSwO++oMICoVB4+HDHvDfW+G2qdpeoHwm05XJ3MS5/Lz1Z+YkzCHTnUmdinUY2nIo/eL60aRKE/3yV3nKNxGIyCgvznHIGPNUHvsTgdhcr+tif/nnNhR4zRhjgE0ishVoBiz04rqFFuRwkJ170rkj+6CiD6qFcouOt+0FX98M057SKatVsTLGsDxpOVM2T+G3bb+RlplGdFg0g5oOYkD8AFrFtNIvf1WggkoEA4FnCjhmBJBXIlgENBaReGAnMBi44ZRjdgC9gLkiUgNoCmwpKOizFeQ8ZYqJI0m+aR84VfPLT0xZHdtZxxeoIktMS+THzT8yZfMUEg8nEh4UTq96vbi0waV0rtWZIEdA9QxXRVTQv5a3jDET8jtARKrktd8Yky0i9wG/YbuPfmKMWS0id3neHw28CIwXkZXYqqTHjTH7C3sT3rJTTORKBIeToNo5vrrcyfo8DzsXw5T7oWars5v2WgW0Y9nH+H3773y/6XsW7lmIIHSq1Ym7295N73q9iQiO8HeIqozKNxEYY0YWdIL8jjHG/AL8csq+0bme7wIuLugaxcWOLHbnXNxTIiihQTLOYNteMLo7fHUzDJsJoRVL5tqqzDLGsDp5Nd9u/JZft/7K4azDxEbGcl/b+7i84eXUqljL3yGqcqCgNgIBBmEbeSdjG3UHAuuA0Xn18CnNThpZnJEGroySqRrKUak2XPMJfHYFTLkPrhnnu4ZqVaalZabxy5ZfmLxxMutS1hHmDKNP/T5c2fhKOtTooPX+qlgVVDX0HlAdCMEmgFDgR2AAtj7/AZ9GV8yCc48s9sWoYm80uAB6PQu/Pwt1Ouj6Beokq/ev5usNX/Pr1l85ln2MZtHNeKrzUwxoMIDIkEh/h6fKqYISQXdjTCsRCQb2ALWMMZki8iWwzPfhFa+g3OsR5IwqjvDD/CldH4CdS+zC97Va28nqVMBKz07n162/8tX6r1idvJrwoHAGxA9gUJNBNK/aXH/9K58rKBFkAxhjskRkUc7UD56GYFf+Hy19gpy5GoszPEMVfDWYLD8icMX78FEvO75g+ByIii3wY6p8SUxL5Kv1X/Hdpu9IzUilYeWGPNHpCS5reJn++lclqqBEsEdEKhpjDhtj+uXsFJGagE/mA/IlWyLwVA1lpNlHfzXYhkbC4C/go4vgqxvhtt8gONw/sagSY4xh4Z6FfLH2C2YnzMYhDi6qdxHXN7te6/6V3xTUa6j/Gd46BFxa/OH4VpBTcBtwuw2OzCN2pz9XEotpDFd9BBMH226lV43RxuNyKsOVwS9bfuGztZ+x8cBGqoRW4Y5Wd3Bt02upWaGmv8NTAc7rUSci0hqIO+Uz3xZ3QL4U7LQLaWS53YRmHrY7Q/zchbNpP7joSZj5kh1f0PV+/8ajitWB9ANMWj+JSesmkZKeQuMqjXnh/BcY0GAAoU6dbkSVDl4lAhH5BGgNrAZyuowaylgiCHLYX9vZLkNohicRhJaCutjuj8KelbYnUfVzoHEff0ekimjHoR18uuZTvt/0PRmuDHrU7cGQ5kPoVLOTVv+oUsfbEsF5xpjmPo2kBDhzEoHbQGYaOEPtQC9/E4ErPoDkLXYxmzt+15HHZdTq/av5eNXH/L79d4IcQVze8HKGNB9Cg6gG/g5NqTPyNhEsEJHmpywqU+bkVA1lu9yQeaR0jewNqQDXT4SPLoQvr7MjjyOi/R2V8kJOA/BHKz/ir91/ERkcye2tbufGc27U5R1VmeBtIpiATQZ7gAzsvEDGGNPaZ5H5QJAzV4kg47B/G4rzEhULg7+E8ZfA10Pgpm8hKMTfUakzMMYwJ3EOH634iBX7VxATHsPD7R9mUJNBVPR325NSheBtIvgEuBlYyYk2gjIn2OFpLHa5IfMwlMa+2rGd4PJ34bvh8MsjcNko7UlUyriNmxk7ZjBmxRjWpayjTsU6PH3e0wxsNFAbgFWZ5G0i2GGMmeLTSEpATonA5TZ2HEFpqhrKrc11kLQO/nwTqjbWnkSlhNu4+X3773zw9wdsOriJuEpxvNT1JQY0GECwoxS0NSl1lrxNBOs800r8iK0aAsAYU6Z6DeU0Fme5jG0j8MeoYm9d9DSkbLHTUETH2zWQlV/klADeX/4+mw5uIr5yPK91f41+cf10oXdVLnibCMKxCSD3lNFlrvtoiPOUqqFKtf0cUT4cDrhyNKQmwjfDYOjPUKe9v6MKKMYY/kj8g3eXv8u6lHXEV47nX93/Rd+4vpoAVLlS0DTU1wPTjDFDSygenwo63mvI01hcGsYQ5Cc43PYkGtsLvhxsu5VWqe/vqALCwt0LeXvZ26xIWkFsZCyvdHuFAfEDNAGocqmgEkF94L+e2UdnAL8CCz1rDJc5wZ42gkyX244jKAs9OypWhxsnw8d94Itr4PZpEJ7nonCqGKxOXs3bS95mwe4F1IiowbNdnmVgo4HaBqDKtYLmGnoNeE1EIoHewG3AaBFZC0wFfjPG7PV9mMUjp2ooO9tV+sYR5KdaUxg80S5oM+kmuPlbCNLeKcUp4VAC7yx7h1+3/UpUaBSPdniUwc0Gay8gFRC8aiMwxqQB33k2RKQ50B/4FOjrs+iKWXCQJxFkZYA7u/SNI8hPXFc7+vib2+G7u+Dqj207giqSA+kH+HDFh3y1/iuCJIhhrYYxtOVQnQZaBZSC2gja5fP2LGPMf4o5Hp/KGVls0j1rEZS1P/ZW18ChnbYnUcUa0O9VHWNwljJcGXy59ks+WvERR7KPcGWjK7mn7T1UjyjhFeuUKgUKKhHk90VvsGsYlxk5bQTu4xPOlZGqodzOvx/S9sD/3odKtexqZ8prxhimb5/Om0veZOfhnXSv052H2z9MoyqN/B2aUn5TUBvBhSUVSEk4XiIoLVNQnw0RuPhlmwymPwMVqkPb6/0dVZmwJnkN/1r4L5buW0rjKo0Z02cMXWp38XdYSvmdt9NQBwN3AzmL684GPjTGZPkoLp/ISQTHVycrS20EueWMMTiaDD/ca3sRNe1X8OcCVEp6CqOWjuLbjd9SJawKz3R5hqsaXaVdQZXy8HZA2QdAMPC+5/XNnn13+CIoX8mpGpIMz+pkpX0cQX6CQu1SlxMug//eAjd/D/X1121u2e5svlr/Fe8te49j2ce4ufnN3NXmLm0IVuoU3iaCjsaYNrlezxSRv30RkC/ldB+VrDJcNZRbaKQdY/BJPzt19dCf7SpniiV7l/DyXy+z8cBGutTqwohOI3RNAKXOwNv+hy4RaZjzQkQaAC7fhOQ7OVVDklmGG4tPVSEGbv7O3stnV8H+Tf6OyK+SjyXz5J9PcuvUWzmceZiRPUfyYZ8PNQkolQ9vSwSPAbNEZAt2LYL6QJmbdiJnHIEjK2fh+nKQCMCuYzDkB1sy+HQg3DbV7gsgbuNm8obJjFw6kmPZxxjWahjDWg8jPCjc36EpVep5O6Bshog0BppiE8E6Y0xGAR8rdXLaCMpdIgCIaWxLBuMvhU8vh6FTIbKGv6MqERsObOD5Bc+zImkFHWt25KnznqJBZS0BKOUtb3sNObEjiOM8n+klIhhj3vRhbMUuZ2EaZ/ZhcIaUv9W/arWGmybDp1fY6Shu+QkqVPV3VD6Tnp3Ohys+ZPyq8USGRPJKt1e4tMGluji8UoXkbdXQj0A6ZXyFModDcDoEZ9bR8lUayC22E9wwCb4Y5EkGU8rlJHWL9iziufnPsSNtB1c0uoJH2j9CVFiUv8NS5VhWVhaJiYmkp6f7O5R8hYWFUbduXYKDvZ8o0dtEUPds1icWkX7A24ATGOuZxO7UY3oCI7HdU/cbYy4o7HUKI9gpBGWX40QAEN8DrvsCJl0Pn19tu5aGVfJ3VMUiLTONN5e8yeQNk4mNjGXsxWPpXKuzv8NSASAxMZHIyEji4uJKbanTGENycjKJiYnEx8d7/Tlvew39KiIXF3zYCZ7qpPewk9M1B673TFaX+5go7NiEy40xLYBBhbnG2Qh2OnC60yEkwteX8q/GvWHQBNj9t52+OmcQXRn2R+IfXPHDFXy78VtubXEr31z+jSYBVWLS09OpWrVqqU0CACJC1apVC11q8TYR/A/4TkSOicghEUkTkUMFfKYTsMkYs8UYkwlMAgaecswNwLfGmB0Axph9hQn+bIQ4HThcGYExjXOzAXDNJ5C42JYMymgySM1I5ck/n+TeGfdSKaQSXwz4gkc6PKI9glSJK81JIMfZxOhtIvgP0AWIMMZUMsZEGmMKqmuoAyTkep3o2ZdbE6CKiMwWkSUiMiSvE4nIcBFZLCKLk5KSvAw5b8FOB05XJgSFFek8ZUbzgbmSQdkrGcxNnMtVP1zFz1t+ZlirYXx16Ve0jGnp77CUKle8TQQbgVWFXJksr7R06ueDgPbAJdheSU+LSJPTPmTMGGNMB2NMh2rVqhUihNMFBwlB7vTAKBHkaHEFXPMxJC6yySC9oMKc/x3NOsrzC57nnhn3UCnUlgLub3c/Ic5y1tNLqUL67rvvEBHWrVtXbOf0NhHsBmaLyBMi8nDOVsBnEoHco5rqArvyOGaqMeaIMWY/8AfQBh8KdjhwugOoRJCjxZW5ksFVkJ7q74jOaNm+ZVw95Wq+2fANQ1sMZdKlk2gR08LfYSlVKkycOJFu3boxadKkYjunt4lgK3bN4hAgMteWn0VAYxGJF5EQYDAw5ZRjfgC6i0iQiEQAnYG13gZ/NoKdDoJMZmCVCHK0uBKunQC7ltsRyMcO+Duik2S5shi1dBS3Tr0VgAn9J/Bwh4d1uUilPA4fPsy8efP4+OOPizUReDuy+PnCntgYky0i9wG/YbuPfmKMWS0id3neH22MWSsiU4EV2PEJY40xqwp7rcKwVUMBWCLIcc5lcN1n8PUQO3Ppzd/b+Yr8bGvqVkbMHcGa5DVc2ehKHu/0OBWCy+g04arce/7H1azZVbxVrM1rV+LZy/Iv+X7//ff069ePJk2aEB0dzdKlS2nXLr+FJL2Tb4lARJ4r6AT5HWOM+cUY08QY09AY87Jn32hjzOhcx7xujGlujGlpjBnpdeRnyZYIAqTX0Jk07Q+DJ8L+jTD+Eji022+hGGP4ZsM3XPfTdew6vIuRPUfyQtcXNAkolYeJEycyePBgAAYPHszEiROL5bwFlQjuKKCbqGCrfJ4rlmhKQLDTQbDJCtwSQY7GveGmb+z01eP62xHIUfVKNITUjFSeX/A807dPp3PNzrzc7WVqVAiM+ZFU2VbQL3dfSE5OZubMmaxatQoRweVyISL8+9//LnK31oLaCD7i5DaBU7eKnmPKjBCng2ATwFVDucV1s7OWHkuxM5cmbSixSy/ft5xrf7yWWTtm8VD7hxhz8RhNAkrlY/LkyQwZMoTt27ezbds2EhISiI+P588//yzyuQtas7jQbQOlXZBTCDEZmghy1O0At/5s1zIY18+WEmqf67PLuY2bcavG8c6yd6hZoSYT+k+gdbVCz16iVMCZOHEiI0aMOGnf1VdfzZdffkn37t2LdG5v5xoqN0IdBiduTQS51Wxl1zD49AoYfxlcPxHii/YPKy8H0g/wxJ9PMG/nPC6ufzHPnv8slULKxxxISvna7NmzT9t3//33F8u5ve0+Wm5UcGTbJ4HcWJyXqg3h9t+gUm07HcXaH4v19Mv3LWfQj4NYuHshT3Z+kjcueEOTgFKlRMAlgnDJtE+0RHC6SrVtyaBWa9u9dPG4Ip/SGMNnaz5j6NShBDuC+WzAZwxuNrhMzNmiVKDwdmGaasAwTixMA4Ax5jbfhOU7EU4tEeQrIto2IP/3VvjpQTi8Dy74J5zFF/eRrCM8O/9Zftv2GxfGXshL3V7SUoBSpZC3bQQ/AHOB3ymDi9bnFuHwhK8lgjMLqQCDv4Qp98PsV+DQTrjkTXB636S0JXULD816iG2HtvFguwe5reVtWgpQqpTy9i87whjzuE8jKSERjiz7REsE+XMGwxXvQ6VaMPc/kLYHBo2zSaIAM3fM5P/+/D9CnaGM6TNG1wxQqpTzto3gJxEZ4NNISki45FQNaYmgQCLQ6xlbGtg03Y5CTtt7xsPdxs17y9/jgVkPEFcpjq8u/UqTgFJlgLeJ4AFsMkj3LErjzcI0pVK42BKBy6HTGXut4+22qihpPYztBftOnxfwcOZhHpj1AKP/Hs3AhgOZ0H8CNSvU9EOwSpVfTqeTtm3b0qZNG9q1a8f8+fOL5bxeJQLPQjQOY0yY57k3C9OUSuGeqqEsh1YNFUrT/jD0F3BlwscXw+ZZx99KOJTATb/cxNzEuYzoNIIXu76oM4Yq5QPh4eEsX76cv//+m1dffZUnnniiWM7rdeufiFwO9PC8nG2M+alYIihhYZ4SQZaEoJVDhVT7XLhjhp2f6POr4ZI3+KtOCx6e/TAiwod9PtSqIKVKyKFDh6hSpUqxnMvb7qOvAR2BLzy7HhCRbsaYEfl8rFQKxSaCDIILXFBB5SEq1o41+OZ2/jvnKV6JiaF+5XjeuehdYivFFvx5pcqDX0fAnpXFe86araD/a/kecuzYMdq2bUt6ejq7d+9m5syZxXJpb0sEA4C2xhg3gIhMAJYBZS4RhHkSQSbaRnC2XCEVeKNROz7PXke3o0d53SFU1PEBSvlcTtUQwIIFCxgyZMjx2UiLojBzDUUBKZ7nlYt0VT8K9YwsTjcBN81SsTiadZTH/niMPxL/4KZzbuIRiSHo13/CR73g+kkQ08jfISrlewX8ci8JXbp0Yf/+/SQlJVG9evUincvbb8NXgWUiMgu7BkEPoHhaKUpYyPGqIS0RFNbeI3u5b+Z9bDywkafPe5prm15r36jWDL6+GT66CAZ9Ao16+zdQpQLAunXrcLlcVK1atcjn8napyokiMhvbTiDA48aYPUW+uh+EGJsItERQOOtT1nPPjHs4nHmYd3u9S7c63U68GdcVhs2CSTfAF4Og17PQ9YGzmpZCKXVmOW0EYOfxmjBhAk6ns8jnzffbUESaGWPWiUjOopiJnsfaIlLbGLO0yBGUsBCTAcAxtyYCb83fNZ+HZz9MheAKfNr/U5pGNz39oCr14bbfYMp98PuzsGsZDHwPQiuWfMBKlVMul29m+Cno2/BhYDjwnzzeM8BFxR6RjwWTRZZxkuEOuIlXz8qUzVN4dt6zxEfF836v9/MfJBZaEa4ZB7XawoznYf8GuO5zO8W1UqrUKmiFsuGep/2NMem53xORMtkNP8idQQbBZGS7/R1KqWaMYezKsYxaNorONTvz1oVvERniRYdbEej2oO0K983tMKYnXPEBnHOpr0NWSp0lb38W5zWOuXjGNpewIJNJOiFkZJfpSVR9ym3cvLbwNUYtG8WA+AF80PsD75JAbo16wZ1/2NLAVzfC9GfBle2bgJVSRVJQG0FNoA4QLiLnYhuKASoBET6OzSeC3JkcIUhLBGeQ6crkyT+fZOq2qdzc/GYe7fAoDjnLarSoerbd4NfHYd5ISFwEV39sZzRVSpUaBbUR9AVuBeoCb+banwb8n49i8imnO5NMo1VDeTmadZQHZz3Igt0LeKj9QwxtMbToawgEhcJlI6FeF7vQzYfd4eqx0KBnMUSslCoOBbURTAAmiMjVxphvSigmn3KabLIIIiNLq4ZyS81I5Z4Z97Bq/ypeOP8Frmx8ZfFeoM11UKuNXQLz0yugx2NwweOFWuxGKeUb3s4++o2IXCIi/xSRZ3I2XwfnCw6TZROBlgiOSzqaxK1Tb2Vt8lrevODN4k8COao3g+GzoO0N8Me/YcJlkLrTN9dSqpzas2cPgwcPpmHDhjRv3pwBAwawYcOGIp3Tq0QgIqOB64B/YNsJBgH1i3RlP3G4s8nEqYnAY+fhnQz5dQg7D+/k/d7v06t+L99eMKSCXfnsyjGw+28Y3RXW/ezbaypVThhjuPLKK+nZsyebN29mzZo1vPLKK+zde+YFo7zhbSvg+caYIcABY8zzQBegTE41Ka5MXBKsvYaAbanbuHXqraRmpjL24rGcV+u8krt4m+tsr6LKsXZE8k8PQ9axkru+UmXQrFmzCA4O5q677jq+r23btnTv3r1I5/W2gjZnDMFREakNJAPxRbqyv7iycEkQmQFeIth4YCPDpg3DYBjXd1zeo4V9LaYR3PE7zHgBFrwL2+fZhuSarUo+FqUK4V8L/8W6lHXFes5m0c14vFP+S8OvWrWK9u3bF+t1wfsSwY8iEgW8DiwFtgETiz2aknC8RBC4iWBdyjpu++02nOL0XxLIERQKfV+Gm76FYwfsxHXzRoE7cP//KFXSCiwRiIgDmGGMOQh8IyI/AWHGmFRfB+cTrkxcEkpGVmB+0azev5rh04cTERzBJxd/UnoWk2nUC+5eAD89ANOfho3TbFtCVD1/R6bUaQr65e4rLVq0YPLkycV+3gJLBJ7FaP6T63WGt0lARPqJyHoR2SQiZ1zERkQ6iohLRK7xKuqicGXhcgQFZBvBiqQVDJs2jMiQSMb1HVd6kkCOClXh2s9g4Puwazm8fz4s+xyM8XdkSpUKF110ERkZGXz00UfH9y1atIg5c+YU6bzeVg1NE5GrpRCji0TECbwH9AeaA9eLSPMzHPcv4Ddvz10krkzcEhJwVUMrklZw5/Q7qRxamXF9x1E3sq6/Q8qbCJx7I9w9z447+OFemDgY0srkrOdKFSsR4bvvvmP69Ok0bNiQFi1a8Nxzz1G7du0indfbxuKHgQpAtoikY7uQGmNMfusTdgI2GWO2eG5gEjAQWHPKcf8AvsGudeB7riyMI7DaCHKSQJWwKnzS95P8ZxAtLarUh1t+hL9G25lM3+sM/f8Nra/VdQ5UQKtduzZff/11sZ7T2wFlkcYYhzEmxBhTyfO6oEVq6wAJuV4nevYdJyJ1gCuB0fmdSESGi8hiEVmclJTkTchn5srE7QgOmJHFK5NWcuf0O4kKjSo7SSCHwwFd7oG7/oSYJvDdcNvV9NBuf0emVLni7YCyGd7sO/WQPPadWtk7ErvaWb7fysaYMcaYDsaYDtWqVSvgsgVwZWKcgVEiWJO85ngSGNdvXNlKArnFNIbbpsLFL8PmmbZ0sPQzbTtQqpjkmwhEJExEooEYEakiItGeLQ4oqFIqkZMHndUFdp1yTAdgkohsA64B3heRKwoRf+G5siAAEsGGAxsYPn04FUMq8nHfj8tuEsjhcML598Hd86FmS7sS2qcDIWWrvyNTAcSUgR8fZxNjQSWCO4ElQDPPY872A7YhOD+LgMYiEi8iIcBgYMopAccbY+KMMXHAZOAeY8z3hb2JQnFlYhzlez2CLQe3MGzaMEKdoXx88cfUrli0hqRSpWpDuOUnuPQt2LkU3u8Cf460CV4pHwoLCyM5OblUJwNjDMnJyYSFFW7dsIJmH30beFtE/mGMeaeQAWWLyH3Y3kBO4BNjzGoRucvzfr7tAj5hDLgycTiDST9aPhNBQloCw6YNQxDGXjy29HURLQ4OB3S4DZr0g18es2skr/wvXPY21O3g7+hUOVW3bl0SExMpcjulj4WFhVG3buF6BXrVa8gY846InA/E5f6MMebTAj73C/DLKfvyTADGmFu9iaVI3C7A4AgO5Uhm+UsEe4/sZdi0YWS4MxjXdxzxlcvmLCBeq1QbBn8Ba3+0CWFsb5sgej0D4VH+jk6VM8HBwcTHl8+/Ka8SgYh8BjQElgM536AGyDcRlDquTAAcQSEczSxfyyampKcwbPowDmYcZOzFY2lcpbG/Qyo551wG8RfA7Fdtd9O1U+Dil6D1ddrVVCkveDuOoAPQ3JTmyjFveBJBUEgoWS5DZrabkKCzXIaxFEnLTOOu6Xex6/AuRvceTcuYlv4OqeSFVYJ+r0KbwXYm0+/uhCUT4JI3oEYLf0enVKnm7bfgKqCMdzsB3LYU4AwOASgXpYIMVwb3z7yfjQc28lbPt+hQM8DryGu1gdunw+XvQNI6GN3drpl87IC/I1Oq1PI2EcQAa0TkNxGZkrP5MjCf8JQIgoNti3pZbyfIdmfzzzn/ZPHexbzU7SW61y3anOTlhsMB7YbAP5ZA+1vgrw/hnfawZLynnUgplZu3VUPP+TKIEpOragjgaEbZLREYY3hhwQvMTJjJiE4juKTBJf4OqfSJiLbdTNsPtaWCHx+ARR9Dv9cgrqu/o1Oq1PB2iok52DUIgj3PF2HXJShbPH3Ngz2JoCyXCN5Z9g7fbfqOO1vfyY3n3OjvcEq3Wq1h6C9w9cdwNAXGD4CvbtbBaEp5eDvFxDDsgK8PPbvqAN/7KCbfyakaKuMlgonrJvLRyo+4uvHV3Nv2Xn+HUzaIQKtr4L5FcOFTsOl3eLcj/Pakth+ogOdtG8G9QFfgEIAxZiNQ3VdB+YwnEYSGlt02gmnbpvHqX6/SM7YnT533FIWYGVwBhETABY/BP5badZMXvAdvt7WP2Rn+jk4pv/A2EWQYYzJzXohIEKdPIFf6eaqGQjyJoKz1Glq6dylPzH2CNtXa8O8e/ybI4W0TjzpNpVow8D24ay7UPhd++z94pwP8/ZUuk6kCjreJYI6I/B8QLiJ9gP8CP/ouLB/xlAjCckoEGWWnRLA1dSv3z7qf2hVr885F7xAeFO7vkMqHmq1gyPdw83d2NPJ3w+HDHrBhms5uqgKGt4lgBJAErMRORPcL8JSvgvKZnKqhcPslWlZKBPuP7efu3+/GKU7e7/0+UWFR/g6p/Gl4EQyfA1eNhcw0+HIQjOsP2xf4OzKlfM7buoVw7KRxH8Hx5SXDgaO+CswnPFVDoSG2RHC4DDQWH8s+xv0z7yf5WDKf9P2E2MhyOIlcaeFwQOtB0HwgLPsU5rwO4/pBw15w4ZNQt72/I1TKJ7wtEczAfvHnCAd+L/5wfMxTInAGhxAe7ORoKW8sdhs3T/75JKv2r+K1Hq/Rqlorf4cUGIJCoOMdcP8y6PMi7FoGYy+CLwfb50qVM94mgjBjzOGcF57nEb4JyYc8iQBHMBVCnRwp5SWCkUtHMn37dB7p8Ai96vXydziBJyQCut4PD66wXU53zIcxPWHi9bBrub+jU6rYeJsIjohIu5wXItIeOOabkHzI5fnidwYTERJUqksE32z4hnGrxnFtk2sZ0nyIv8MJbKGRtsvpgyttFdH2eTDmAvjiWkhc7O/olCoyb9sIHgT+KyI5S03WAq7zSUS+5Jl0DoeTiJDSWyJYtGcRL/3vJc6vfT4jOo/QsQKlRVhluOCf0PlOWDjGjj0Y2wsaXAg9HoX6XXXaa1UmebswzSIRaQY0xS5Kv84YU/bWBnR7QnYEUyG0dJYIEg4l8NDsh4itFMsbF7xBsCPY3yGpU4VVhh6PQee77NxFC96D8ZdA3U7Q/RFofLFteFaqjCjMv9aOQGvgXOB6ESl79RU569o6g22JoJR1H03LTOO+mfcB8O5F7xIZEunniFS+QiOh24O2DWHAG5C2ByZeBx+cD8sn6jrKqszwdq6hz4A3gG7YhNARu1hN2XK8aiiICiFBHC1FA8pcbhf//OOf7Di0g7d6vkW9SvX8HZLyVnA4dBoG9y+FK8eAOOD7u+zUFfPfgfRD/o5QqXwF1gpl7lyNxaGlq0Tw9tK3+XPnnzx93tN0rNnR3+Gos+EMtvMXtb4WNk6H+aNg2lMw5992fYTOd0GUjgNRpU9grVCWU1TPKRGUkjaCn7b8xLjV47iu6XVc2/Raf4ejikoEmlwMt/4Ew2dD4z7wvw/g7Tbw9S2QsFCnr1ClirclgpwVyhYCx6doNMZc7pOofOV41ZCnRFAKeg2tSV7Dc/Ofo0ONDjze6XF/h6OKW+1z4ZpPoPfztqfRkgmw5nuo3c6WEFpcAUGh/o5SBbgAW6Esp0TgpEJIEBnZbrJdboKc/unhkZKewoOzHiQ6LJr/9PyP9hAqz6Ji4eIX4YLH4e+JdvnM74bDtCeh/a12FbXKdfwdpQpQhVmhbB0Q6dnWevaVLe4scASDCBEhTgCOZvmneijbnc2jcx4lJT2FkReOJDos2i9xqBIWWtE2LN+7EG76Fuq0hz/egJGt4KubYMtsnQZblTivSgQici3wOjAbO47gHRF5zBgz2YexFT93Nnjm8K8Qah+PZrioFFbyv8TfWvIWi/Ys4uVuL9O8avMSv77yM4cDGvWy24FtdjzCss9g7Y8Q3RA6DIU2N0CFqv6OVAUAb+tEngQ6GmNuMcYMAToBT/suLB9xZdueHXC8ROCPnkNTt07l0zWfcn2z67m8YdlqZlE+UCXOVhs9vM52P61QzfY2erMZTL4Ntv6hjcvKp7xtI3AYY/blep1M4QajlQ7urBMlgpATJYKStOnAJp6Z/wxtq7XlsQ6Plei1VSkXHGa7n7a5DvaugaUTbHvCqm+gSjycexO0vdGurqZUMfL2y3yqiPwmIreKyK3Az8CvvgvLR1xZJ0oEobZEkJZRcqM/j2Qd4aHZDxERFGEbh53aOKzOoEZz6P8veGS9LSVUrgszX4S3msMXg2D197rGsio23s419JiIXIUdWSzAGGPMdz6NzBfc2baxGIgKDwHg0LGSSQTGGJ6Z9wwJaQl8dPFHVI+oXiLXVWVccPiJUkLyZlj2Ofw9Cf57C4RXgZbXQJvroU47nfBOnbV8E4GINAJqGGPmGWO+Bb717O8hIg2NMZtLIshi484Ghy0JRFewiSDlSMkkgi/Xfcm07dN4qP1DOnJYnZ2qDaH3s3DRU7BlFiz/0jYwL/oIqja2yaLVINvmoFQhFFQ1NBJIy2P/Uc97+RKRfiKyXkQ2iciIPN6/UURWeLb5ItLGm6DPWq6qoagI+3jgaKZPLwmwImkFbyx+g56xPbm1xa0+v54q5xxOaNTbDlR7dANc/g5UrA4zX7Kjlz/uC4vGwpFkf0eqyoiCqobijDErTt1pjFksInH5fdCzrvF7QB8gEVgkIlOMMWtyHbYVuMAYc0BE+gNjgM6FuYFCyRlHAIQFO6kQ4iTliG8TQWpGKo/OeZQaETV4qetLOKTstbGrUiyssp3HqN0QOLgDVv4XVnwNPz8Cvz5u10poeTU0G2CPVSoPBSWCsHzeC8/nPbBdTDcZY7YAiMgkYCBwPBEYY+bnOv5/QN0Czlk0rmxwnrjlKhVCOODDRGCM4ck/nyTpWBKf9f+MyqH6h6h8KKqeXQ+h28Owd7VNCqu+sTOhOkNtKaLlVdCkr51CWymPghLBIhEZZoz5KPdOEbkdWFLAZ+sACbleJ5L/r/3bOUNPJBEZDgwHqFevCNMz5yoRgG0nSPFh1dCnaz5lTuIcRnQaQcuYlj67jlInEYGaLe3W+zlIXASrvrVzHK3/+URSaD4QmvbTkoIqMBE8CHwnIjdy4ou/AxACXFnAZ/PqwpDnqBgRuRCbCLrl9b4xZgy22ogOHTqc/ciaXCOLAapE+K5EsDJpJSOXjOSi2Iu4odkNPrmGUgUSgdhOduv7CiT8ZRPCmh9sUnAEQ4OecM5l0HQAVKzm74iVH+SbCIwxe4HzPV/UOT9pfzbGzPTi3IlA7snX6wK7Tj1IRFoDY4H+xhjftm7lGlkMULVCCFv2Hy72y6RlpvHYH49RLaIaL3R9QdccVqWDwwH1u9it76uwczGsnQJrpsCP98NPD0LsedDsEmja3/ZSUgHB23EEs4BZhTz3IqCxiMQDO4HBwEk/jUWkHrZL6s3GmA2FPH/hubMg6ESzR5UKIaQcLt4SgTGGFxa8wJ4jexjfb7y2C6jSyeE4UVLo8yLsXQVrf4J1P9kZUac9CdWa2YTQpD/U7XC867Uqf7ydYqLQjDHZInIf8BvgBD4xxqwWkbs8748GngGqAu97fjVnG2N8twTmKVVD0RVCOJLpIj3LRVhw8fwj/37T90zdNpUH2j1A2+pti+WcSvmUCNRsZbcLn4AD22H9r7bqaP478OdbEFEVGvWxC+407AXhUf6OWhUjnyUCAGPML8Avp+wbnev5HcAdvozhJLnGEYBtIwA4eDSLmpWLngi2pm7l1YWv0rlWZ25reVuRz6eUX1SpD+fdZbdjB2HzDFg/FTb+BismgTghtjM07m2TQ81WOqq5jPNpIih13K5TSgQ2KaQcyaRm5fx6yhYs05XJ4388TqgzlFe6vaLjBVT5EB5lxyG0vNr+/SQutglh43SY8YLdKtaAhhfZrcGF2uBcBgVYIji91xAUz+jiUUtHsTZlLaMuHKXzCKnyyeGEep3t1usZSNtrSwubfocNv9mZUsGWEHKSQr3z7HxJqlQL6ERwYr6hoiWCBbsWMGHNBAY1GcSF9S4s0rmUKjMia0DbG+zmdsHu5bB5JmyeDQveh3lv2zELsZ2gwQUQf4Fdw1ln3S11AjoRVKlQ9BJBakYqT/35FHGV4niso64voAKUw2mX3azTHno8BhmHYft82DoHtsyx8yDxEoRUhHpdIL47xHWDmm1OGu2v/COw/g+c0kYQFX6ijeBsGGN4fsHzpKSn8E6vdwgP0iKwUoBdm7nJxXYDOLIftv1pV1vbNhemT7f7Qyrahue4rlC/qy0xBIX6L+4AFWCJIPukvtBBTgeVw4PPenTxT1t+Yvr26TzQ7gFdd1ip/FSIgRZX2A1s+8L2eTY5bJ9vG53BjvOp096WGup1gdiOOgVGCQiwRJB1UokAbDvB/rNIBLsP7+aVv16hXfV2DG0xtLgiVCowRNawE+C1vMq+PrIfdiyA7Qtgx3w7dsG8AQhUb35i8FtsZ4huoN1Vi1mAJYLs0xJBrcph7Dp4rHCnMW6emvcUbuPm5W4v49QRl0oVTYUYO9/ROZfZ1xmHYecS2PE/SPifnUV1yTj7Xng01O3o2dpD7XY6wK2IAiwRuE7rsVAvOoLf1+4r1Gm+XPslC/cs5Pnzn6dupG9nzlYqIIVWtD2NGlxgX7vdsH+9nTQvcREkLLLjGXLENLFVSrXb2WU7a7SE4KKNDQokAZYIsk+bLyU2OoL9hzM4mplNREjB/zm2pG5h5NKR9KjbgysbFTQBq1KqWDgcUP0cu7W/1e47dhB2LYXEJXYCvU0zToxlcATZKqXa50LttlCrLdRooQ3RZxCAieDkW65bxfb0STxwjCY18l+sI9udzVN/PkVYUBjPdXlOZxVVyp/Co06MaAYwBg7thJ1LYdcyu635AZZOsO87gqF6M9tltVZrqNnartmgi/QEUCIwJs9EEBsdAUBCytECE8H41eNZuX8lr/d4nWoROoxeqVJFBCrXtVvzy+0+Y+Dgdti13A542/03bPgVln9+4nPRDWxVUs1WnseWUDk2oBqkAygRuO3jqYmgyolEkJ8NBzbw3vL36FO/D/3i+/kkRKVUMROBKnF2y+m6agwc2gV7Vnq2FZ5puKec+FxoJVu1VKO557GFrZYKr+KHm/C9wEkEriz7eEobQUzFEMKDnexIOXPPoSx3Fk/9+RSVQirx1HlP+TJKpZSviUDlOnZrmutHXUYa7F0D+1bbNZ/3roaV30DGJyeOqVjzRFtFtaZQzfNYxnstBU4icGfbx1NKBCJCbHQ4CQfOXCIYt2oca1PW8mbPN4kOi/ZllEopfwmNPDGpXo6c0sO+NZ5tnX1cMh6ycn1nVKxhey5Va2ofc7ZKtctEFVMAJoLTJ7yKrRJxxqqhjQc28sHfH9Avrh996vfxZYRKqdImd+mhca6/f7cbUnfYxLB/PSRtsI8r/gsZqSeOC65gl/yMaQxVG3m2hhDdsFSVIgIoEbjsYx6Dv2KjI/hrawrGmJN6AmW7s3l63tNEBkfyROcnSipSpVRp53CcaHvIXb1kDBzeZ5PC/o12S95oxz6s+hYwJ46NqGoTQtWGtsE6ugFEx0OVeIgo2ZqHAEoEeVcNgR1Udjgjm6S0DKpXOjEIZcLqCaxOXs3rF7yuVUJKqYKJ2OkzImtAfI+T38tKhwPbIHkTpGyG5M2QssVOxJcz/iFHWGWbEKrEeZJDHETVt20TkTWLPewATASnlwha1bWTWv2dmEqf5jYRbE3dyvvL36d3vd70rd+3xMJUSpVTwWF2HEP1Zqe/l3nUJokDW21ySNlqX+9ZAet+OvH9df4/4OKXij20wEkEJqdq6PRbblm7Mk6HsDzhAH2a18Bt3Dw7/1nCgsJ48rwndeCYUsq3QiJsV9Uaecxi7HbZgXIHtttGaR8InESQT9VQeIiTZjUjWZ5wEICJ6yaybN8yXu72MjHhMSUYpFJKncLhhKh6dvPVJXx25tLGfeYSAUDb2ChWJKSyIzWBt5e+Tbc63biswWUlGKBSSvlHACWCM7cRgE0EaRlZjPjjaRzi4Nkuz2qVkFIqIAReIpC8E8G59aIIjlrIypQlPNz+YWpWKP6WeaWUKo0CKBHkXzUUHp5GWI1fqMQ5DGoyqAQDU0op/9JEQM4i9M8R5IB9Wy7j4NGsEg5OKaX8J4ASwZnbCCZvnMyC3Qu4uem9ZGVE8/PK3SUcnFJK+U/AJ4Idh3bw+qLXOa/WeTzQaQiNq1dk8pJEjDF5nEQppcqfwEkEeQwoy3Jn8cTcJwhyBPFi1xdxOpwM6VKf5QkHmbpqj58CVUqpkhU4iSCPAWWjlo5ixf4VPNPlmeO9hK7vVI9mNSN58ac1HM3M9kekSilVogIoEXhKBJ7uo9O2TWP86vFc1/Q6+sWdmD0wyOnghYEt2X0oneGfLiE9y+WPaJVSqsT4NBGISD8RWS8im0RkRB7vi4iM8ry/QkTa+SyYXNNQL9+3nCfmPkHbam15rONjpx3aKT6a169pw7zN+7nuwwWs2pl62jFKKVVe+GyuIRFxAu8BfYBEYJGITDHGrMl1WH+gsWfrDHzgeSx+nqqhGfuX8cSqD6lZoSajLhpFqDM0z8OvaV+XiBAnT3+/ikvf+ZOOcVXo2bQ6repUpl50BNUiQ4kIceroY6VUmefLSec6AZuMMVsARGQSMBDInQgGAp8a20XnfyISJSK1jDHF3n9z3oG1vFi3NjuXj+Sc6HN4r9d7VAnLfyHqAa1q0bVhDF8s3M4Py3bx+m/rT3o/yCGEhzgJDXIS7BScDrs5RDieHoTjz/NLGppOlFIFua5jLHd0b1Ds5/VlIqgDJOR6ncjpv/bzOqYOcFIiEJHhwHCAevXObga+CpXq0Coshltb3MjVrW8n2Hn6kpV5qRwRzD09G3FPz0YcPJrJuj1pJKQcJflIJmnpWRzJcJHpcpOV7cblNriMwe3peWqMObEeUT69UU1+byqllEdMxbxrMIrKl4kgrx+5p37jeXMMxpgxwBiADh06nNW3ZtuWN9C25Q1n89HjoiJCOK9BVc5rULVI51FKqdLEl43FiUBsrtd1gV1ncYxSSikf8mUiWAQ0FpF4EQkBBgNTTjlmCjDE03voPCDVF+0DSimlzsxnVUPGmGwRuQ/4DXACnxhjVovIXZ73RwO/AAOATcBRYKiv4lFKKZU3ny5VaYz5Bftln3vf6FzPDXCvL2NQSimVv8AZWayUUipPmgiUUirAaSJQSqkAp4lAKaUCnJS1BVhEJAnYfpYfjwH2F2M4ZYHec2DQew4MRbnn+saYanm9UeYSQVGIyGJjTAd/x1GS9J4Dg95zYPDVPWvVkFJKBThNBEopFeACLRGM8XcAfqD3HBj0ngODT+45oNoIlFJKnS7QSgRKKaVOoYlAKaUCXMAkAhHpJyLrRWSTiIzwdzy+ICKxIjJLRNaKyGoRecCzP1pEpovIRs9j/mt0ljEi4hSRZSLyk+d1eb/fKBGZLCLrPP+vuwTAPT/k+Te9SkQmikhYebtnEflERPaJyKpc+854jyLyhOf7bL2I9C3KtQMiEYiIE3gP6A80B64Xkeb+jconsoFHjDHnAOcB93rucwQwwxjTGJjheV2ePACszfW6vN/v28BUY0wzoA323svtPYtIHeB+oIMxpiV2WvvBlL97Hg/0O2Vfnvfo+bseDLTwfOZ9z/fcWQmIRAB0AjYZY7YYYzKBScBAP8dU7Iwxu40xSz3P07BfEHWw9zrBc9gE4Aq/BOgDIlIXuAQYm2t3eb7fSkAP4GMAY0ymMeYg5fiePYKAcBEJAiKwKxmWq3s2xvwBpJyy+0z3OBCYZIzJMMZsxa7p0ulsrx0oiaAOkJDrdaJnX7klInHAucBfQI2cld88j9X9GFpxGwn8E3Dn2lee77cBkASM81SHjRWRCpTjezbG7ATeAHYAu7ErGU6jHN9zLme6x2L9TguURCB57Cu3/WZFpCLwDfCgMeaQv+PxFRG5FNhnjFni71hKUBDQDvjAGHMucISyXyWSL0+9+EAgHqgNVBCRm/wbld8V63daoCSCRCA21+u62KJluSMiwdgk8IUx5lvP7r0iUsvzfi1gn7/iK2ZdgctFZBu2uu8iEfmc8nu/YP8tJxpj/vK8noxNDOX5nnsDW40xScaYLOBb4HzK9z3nONM9Fut3WqAkgkVAYxGJF5EQbCPLFD/HVOxERLB1x2uNMW/memsKcIvn+S3ADyUdmy8YY54wxtQ1xsRh/5/ONMbcRDm9XwBjzB4gQUSaenb1AtZQju8ZWyV0nohEeP6N98K2f5Xne85xpnucAgwWkVARiQcaAwvP+irGmIDYgAHABmAz8KS/4/HRPXbDFg9XAMs92wCgKrbHwUbPY7S/Y/XBvfcEfvI8L9f3C7QFFnv+P38PVAmAe34eWAesAj4DQsvbPQMTsW0gWdhf/Lfnd4/Ak57vs/VA/6JcW6eYUEqpABcoVUNKKaXOQBOBUkoFOE0ESikV4DQRKKVUgNNEoJRSAU4TgVJKBThNBCpgiUhVEVnu2faIyE7P88Mi8r4PrjdeRLaKyF35HNNdRNbknopYKV/TcQRKASLyHHDYGPOGD68xHjvobXIBx8V5jmvpq1iUyk1LBEqdQkR65lrk5jkRmSAi00Rkm4hcJSL/FpGVIjLVM7cTItJeROaIyBIR+S1nfpgCrjPIs9DK3yLyh6/vS6kz0USgVMEaYtc8GAh8DswyxrQCjgGXeJLBO8A1xpj2wCfAy16c9xmgrzGmDXC5TyJXygtB/g5AqTLgV2NMloisxK6ONdWzfyUQBzQFWgLT7ZxoOLFzxhRkHjBeRL7GzqiplF9oIlCqYBkAxhi3iGSZEw1rbuzfkACrjTFdCnNSY8xdItIZW9pYLiJtjTHJxRm4Ut7QqiGlim49UE1EuoBdE0JEWhT0IRFpaIz5yxjzDLCfk+eXV6rEaIlAqSIyxmSKyDXAKBGpjP27GgmsLuCjr4tIY2yJYgbwt08DVeoMtPuoUiVEu4+q0kqrhpQqOanAiwUNKAN+xFYVKVUitESglFIBTksESikV4DQRKKVUgNNEoJRSAU4TgVJKBbj/ByWHgC+8R0/cAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sp=0 # cell number 0, camp only saved cell No 1 \n",
    "plt.figure()\n",
    "plt.plot(solTchem[:,sp,t_indx],solTchem[:,sp,A_indx], label=\"A\")\n",
    "plt.plot(solTchem[:,sp,t_indx],solTchem[:,sp,B_indx], label=\"B\")\n",
    "plt.plot(solTchem[:,sp,t_indx],solTchem[:,sp,C_indx], label=\"C\")\n",
    "plt.xlabel('Time [s]')\n",
    "plt.ylabel('Concentration [mol/m3]')\n",
    "plt.legend(loc='best')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7fa8501e9970>"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAmO0lEQVR4nO3deZxU1Zn/8c/TVdXdIA0YFlFQQeJoEGRr3GZ0XMa4ixhx0EkmbuPy0wlOll90ohnioHFJnLiG4J7IQEYUt5dxHTT+YlRQkFUUFbVVFnFh0e6mm+f3x73VFE13123oW9Xd9/t+vfrVd616Ton19DnnnnPM3RERkeQqKXYAIiJSXEoEIiIJp0QgIpJwSgQiIgmnRCAiknDpYgfQWr179/aBAwcWOwwRkQ7ltdde+9Td+zR1rsMlgoEDBzJ37txihyEi0qGY2fvNnVPTkIhIwikRiIgknBKBiEjCdbg+AhGRrE2bNlFVVUV1dXWxQ2k3ysvLGTBgAJlMJvI9SgQi0mFVVVVRUVHBwIEDMbNih1N07s7atWupqqpi0KBBke+LrWnIzO42s9VmtqiZ82ZmN5vZcjNbYGaj4ooFgBduhEk94C83x/o2IlI41dXV9OrVS0kgZGb06tWr1TWkOPsI7gWObeH8ccDe4c/5wG9jjAVm/yL4/cyVsb6NiBSWksDWtufziK1pyN3/bGYDW7hkLPB7D+bBftnMeprZru7+SZsGMqlH88cmfdmmbyUi0hEV86mh/sCHOftV4bFtmNn5ZjbXzOauWbOmde9S0kyua+64iEgrrFy5kgkTJjB48GCGDBnC8ccfz1tvvQXAf/3Xf1FeXs6XX275o/P555/HzLjrrrsajs2bNw8z41e/+hUAZ511FoMGDWLEiBGMGjWKv/71r7GWoZiJoKn6S5Or5Lj7VHevdPfKPn2aHCHdvH9bvNULN7zBvy1p3euIiDTi7owbN47DDz+cd955hyVLlnDNNdewatUqAKZPn86YMWOYNWvWVvcNGzaMP/7xjw37M2bMYPjw4Vtdc8MNNzB//nyuvfZaLrjggljLUcxEUAXsnrM/APi4zd+lol/DZr3n5J6KXdr8rUQkWWbPnk0mk+HCCy9sODZixAgOPfRQ3nnnHTZs2MDkyZOZPn36VvftscceVFdXs2rVKtydJ598kuOOO67J9zjssMNYvnx5rOUoZvvIo8AlZjYDOBD4ss37B7L2PZGVm3twzsKh3LP/YvrZF7G8jYgUzy8eW8ySj9e16WsO2a07/3HSfs2eX7RoEaNHj27y3PTp0znjjDM49NBDWbZsGatXr6Zv374N50877TQeeOABRo4cyahRoygrK2vydR577DGGDRu2YwXJI87HR6cDfwX2MbMqMzvXzC40s2zqfAJ4F1gO3AH8n7hiYcI0qg6ZzKfenS6fvwUn3BjbW4mIQNDcM2HCBEpKSjj11FN54IEHtjp/+umn88ADDzQkjMZ+8pOfMGLECKZOnbpVf0Ic4nxqaNuSbX3egYvjev/GMqkSfpCeRffVc+CF6+BEJQORzqSlv9zjst9++zFz5sxtji9YsIC3336bo48+GoDa2lr22msvLr54y1dev379yGQyPPPMM9x000289NJLW73GDTfcwGmnnRZvAULJmGtocl9G3LUn30s/i+Ew967gEdLJffPfKyLSjCOPPJKamhruuOOOhmNz5sxh4sSJTJo0iRUrVrBixQo+/vhjPvroI95/f+uZoK+66iquu+46UqlUoUPfSjISwcQFrNv7FL720mA/3QWGjYeJC4sbl4h0aGbGrFmzeOaZZxg8eDD77bcfkyZN4vnnn2fcuHFbXTtu3DhmzJix1bFDDjmEU045pYARN82CFpqOo7Ky0rdnYZp1My+h28L78VQpKd8Eo89W85BIB7d06VK+9a1vFTuMdqepz8XMXnP3yqauT8yoqszXn3J//VHsctiFHPP1n2DDqmKHJCLSLiQmEawfey8/v+Y5Ju/0N3DU0cUOR0Sk3UhGHwHBU0MAm+o3FzkSEZH2JTmJIK1EICLSlOQkglQwvcSm+o7VOS4iErfkJIKSEvrwOWPnnQfr1VEsIpKVmERQUmJMzMxit3Xzg5HFIiI7aO3atYwYMYIRI0bQr18/+vfv37B//fXXs++++zJ06FCGDx/O73//ewAGDhzIp59+WuTIt5aMp4Ym94W6Gr6bHbw3967gJ10GV6wuamgiUmDrV8LMs+G0e3d4FuJevXoxf/58ACZNmkS3bt348Y9/zJQpU5g1axavvvoq3bt358svv+Thhx/e4dDjkowawcQFMHQ8X6ORxSKJ98L18MHLsbYMXHPNNdx+++10794dgB49evD973+/4fwtt9zCqFGjGDZsGG+++SYAGzdu5JxzzmHMmDGMHDmSRx55BIB7772XU045hZNOOolBgwZx6623cuONNzJy5EgOOuggPvvssx2ONxmJoKIflFVQxiY2WSnU10BZd61JIJIkk/sGc4zNvQt8c2xzjq1fv57169czePDgZq/p3bs3r7/+OhdddFHDqmRXX301Rx55JHPmzGH27Nn85Cc/YePGjUAw3fV///d/8+qrr/Kzn/2Mrl27Mm/ePA4++OCGJqcdkYxEALBxNQ+VHMOte00JppfQyGKRZAlbBkh3CfZjahlw97wLyJ966qkAjB49mhUrVgDw9NNPc+211zJixAgOP/xwqqur+eCDDwA44ogjqKiooE+fPvTo0YOTTjoJCFY6y96/I5LRRwAwYRo3Xz+b0aU7w4nfKXY0IlJoYcsA9TWQLo+tZaB79+7stNNOvPvuu+y1115NXpNdhCaVSlFXVwcECeTBBx9kn3322eraV155ZatFa0pKShr2S0pKGu7fEcmpEQDplFGrAWUiybVxddAicN6zsbYMXH755Vx88cWsWxesmLZu3TqmTp3a4j3HHHMMt9xyC9mJQOfNmxdLbE1JTo0AKE2VsKlOiUAksSZM27Id4+zDF110ERs2bGDMmDFkMhkymQw/+tGPWrznyiuv5NJLL2X//ffH3Rk4cCCPP/54bDHmSsw01AAn3fL/6N2tlHvOPqCNoxKRYtA01E1r7TTUiWoayqRMU0yIiDSSsERQoknnREQaSVQiKE0rEYh0Nh2teTtu2/N5JCoRBDUC/aMR6SzKy8tZu3atkkHI3Vm7di3l5eWtui9RTw0FfQSqEYh0FgMGDKCqqoo1a9YUO5R2o7y8nAEDBrTqnoQlghKNIxDpRDKZDIMGDSp2GB1ei4nAzE6N8BrV7v5EG8UTq1J1FouIbCNfjeAO4BGgpYkzDgM6RCLIpErYVKe2RBGRXPkSwZ/c/ZyWLjCz+9swnlil1UcgIrKNFp8acvfv5nuBKNe0Fztv/ozf1V+ppSpFRHLkfXzUzPqZWb9wu4+ZnWpm+8UfWts7YtW9jPI3tVSliEiOfJ3FFwCXBZt2HXAWsBj4pZld7+53xR9iGwiXqhwNQW+HlqoUEWmQr4/gEmA/oAvwPvBNd19pZjsDs4GOkQgmLoCnrmDTkkfJbK7B012wb50I37662JGJiBRdvqahTe7+lbuvBd5x95UA7v45kPfxGzM71syWmdlyM7usifM9zOwxM3vDzBab2dnbVYp8wgUp0ptrqfaMlqoUEcmRLxFsNrNMuH1C9qCZlee718xSwG3AccAQ4AwzG9LosouBJe4+HDgc+LWZlUYPvxU2rmZJ//GMq72KTSPP0lKVIiKhfE1DpxL+5e/uVTnHewEtr7IABwDL3f1dADObAYwFluRc40CFBQt8dgM+A3Z83bWmTJjGnL+8x9J3lrDxqPMo3SmefCMi0tHke3z0A3dv+GI2s+5m9g3ga+D1PK/dH/gwZ78qPJbrVuBbwMfAQmCiu2/zoL+ZnW9mc81s7o7MKZJJB8XVWAIRkS0izT5qZheY2SpgAfBa+JNvmbCmRiM37lc4BpgP7AaMAG41s+7b3OQ+1d0r3b2yT58+UUJuUiYVFFfzDYmIbBF10rkfA/u5+6eteO0qYPec/QEEf/nnOhu41oM5ZJeb2XvAvsCrrXifyEpT2RqBppkQEcmKuh7BO8BXrXztOcDeZjYo7ACeADza6JoPgKMAzGwXYB/g3Va+T2TpVFBJUdOQiMgWUWsElwMvmdkrQE32oLv/oLkb3L3OzC4BngJSwN3uvtjMLgzPTwH+E7jXzBYSNCX9tJW1jlZpaBqqUyIQEcmKmgh+B/wvQYdu5G/RcHrqJxodm5Kz/THw7aivt6O2NA0pEYiIZEVNBHXu/sNYIymAjPoIRES2EbWPYHb4COeuZvaN7E+skcUgoz4CEZFtRK0RnBn+vjznmAN7tW048dI4AhGRbeWbfXRXd//E3TvFoqB6fFREZFv5agR3hzONPg88Cfy/3JHGHU1GncUiIttoMRG4+3HhBHOHA+OAX5nZBwRJ4Ul3/yD+ENuO+ghERLaVt4/A3asJv/gBzGwQwYyit5pZP3c/IN4Q247GEYiIbCtqZ3EDd38PuB24PbYpo2NSmlYfgYhIY/k6i9ezZaK47CRyHm67u28zQVx7pj4CEZFt5esjqChUIIWguYZERLYVuWnIzIYDh4a7f3b3BfGEFJ9STUMtIrKNqOsRTASmAX3Dn2lm9q9xBhaHzFer+WPpVWS+2v7FbUREOpuoNYJzgQPdfSOAmV0H/BW4Ja7A4pB68QbG2DJK378DOLjY4YiItAtRE4EB9Tn79TS9Aln7NLkv1AWzZ5cYjFz1IEx6ENJlcMXqIgcnIlJcUSeduwd4xcwmmdkk4GXgrtiiamsTF8DQ8ZDuAkCtlcGw8TBxYZEDExEpvkg1Ane/0cyeB/6OoCZwtrvPizOwNlXRD8oqoL6GGjJkvBbKukPFLsWOTESk6FozoOw9oC68x8xslLu/Hk9YMdi4Gkafzbnzh3Bx9xc5eMOqYkckItIuREoEZvafwFkEaxdnB5g5cGQ8YcVgwjQAViz+Xx7oV8nBp48objwiIu1E1BrB6cBgd6+NM5hCKE2VaIoJEZEcUTuLFwE9Y4yjYDKpEjZp0jkRkQZRawS/BOaZ2SKgJnvQ3U+OJaoYpVOmKSZERHJETQT3AdcBC4EO/S2aSZVoigkRkRxRE8Gn7n5zrJEUSNBHoEQgIpIVNRG8Zma/BB5l66ahjvP4aCiTNqo3KRGIiGRFTQQjw98H5RzrWI+PhjKpEjZUd9hll0VE2lzUkcVHxB1IoQR9BHp8VEQkq8XHR83sxHwvEOWa9kR9BCIiW8tXI7jBzD6i5ZlGrwEeb7uQ4pXR46MiIlvJlwhWATfmuebtNoqlIDSgTERka/nWLD68QHEUTCatPgIRkVxRp5jYLmZ2rJktM7PlZnZZM9ccbmbzzWyxmb0QZzygPgIRkcZaMw11q5hZCrgNOBqoAuaY2aPuviTnmp7A7cCx7v6BmfWNK56sdIn6CEREcm13jcDMMnkuOQBY7u7vhrOWzgDGNrrmTOAhd/8AwN1jXzcyk1aNQEQkV6sSgQWONLM7Cf7Kb0l/4MOc/arwWK6/AXY2s+fN7DUz++dm3vd8M5trZnPXrFnTmpC3kQmnoXZXP4GICERMBGZ2oJndBLxPMM3Ei8C++W5r4ljjb980MBo4ATgGuNLM/mabm9ynunulu1f26dMnSsjNKk0FYWlNAhGRQL4BZVeb2dsEYwUWEkw1scbd73P3z/O8dhWwe87+AODjJq550t03uvunwJ+B4a0pQGtlUkGR1TwkIhLIVyM4n2AswW+B+919Ldv+Vd+cOcDeZjbIzEqBCQS1iVyPAIeaWdrMugIHAksjR78dsomgTjUCEREg/1ND/YBvA2cAvzGz2UAXM0u7e4szt7l7nZldAjwFpIC73X2xmV0Ynp/i7kvN7ElgAcE6B3e6+6IdLFOLMukgEWhNAhGRQL4BZfXAn4A/mVk5cCLQFfjIzJ5z9zPz3P8E8ESjY1Ma7d8A3LAdsW+XLX0ESgQiItCKcQTuXg3MBGaaWXdgXGxRxUh9BCIiW9uuAWXuvo5g+coOR4lARGRrsU4x0R5lE0FtnTqLRUQg/+Oj48PfgwoTTvy6bVrDH0uvwtevLHYoIiLtQr4aweXh7wfjDqRQ9lp8O2NsGX1ev6nYoYiItAv5+gjWho+MDjKzxmMAcPeT4wkrBpP7Ql0NuwEY9F12P0y6H9JlcEXsUxyJiLRb+RLBCcAo4A/Ar+MPJ0YTF8BTV1C/9DFS9dXUp8pJDTkJvn11sSMTESmqfOMIaoGXzewQd19jZhXBYd9QmPDaUEU/KKugpL6Gas9QVl8DZd2hYpdiRyYiUlRRnxraxczmAYuAJeFMoUNjjCseG1fz+ZDvMa72Kj7YawJsWFXsiEREii7qOIKpwA/dfTYEq4qFxw6JJ6yYTJjGZ6s3sPT1F3hj+CnsOXy3YkckIlJ0UWsEO2WTAIC7Pw/sFEtEMSvNDijTAvYiIkD0GsG7ZnYlQacxwHeB9+IJKV6ZtOYaEhHJFbVGcA7QB3go/OkNnB1XUHHSFBMiIluLVCMIF6H5QcyxFETDFBNaj0BEBEjgXEOlqhGIiGwlcYkgnV2PQJ3FIiJAEhNBiTqLRURyReojMLM+wL8AA3Pvcfdz4gkrPmZGaapEfQQiIqGoj48+ArwIPAvUxxdOYWRSphqBiEgoaiLo6u4/jTWSAsqkS5QIRERCUfsIHjez42ONpIAyqRI2qWlIRASInggmEiSDajNbH/6sizOwOJWmVCMQEcmKOqCsIu5ACkl9BCIiW0TtI8DMTgYOC3efd/fH4wkpfhnVCEREGkRqGjKzawmah5aEPxPDYx1SJlVCbZ36CEREIHqN4HhghLtvBjCz+4B5wGVxBRYnPTUkIrJFa0YW98zZ7tHGcRRUpkR9BCIiWVFrBL8E5pnZbMAI+goujy2qmKmPQERki6hPDU03s+eBMQSJ4KfuvjLOwOKUSZfw1dcdfoC0iEibaLFpyMz2DX+PAnYFqoAPgd3CYx1Saco0+6iISChfjeCHwPnAr5s458CRbR5RAahpSERkixYTgbufH24e5+7VuefMrDzfi5vZscBNQAq4092bfOTUzMYALwP/6O4zowS+IzKpEuo26/FRERGI/tTQSxGPNTCzFHAbcBwwBDjDzIY0c911wFMRY9lhwTgC1QhERCBPjcDM+gH9gS5mNpKgoxigO9A1z2sfACx393fD15oBjCUYkJbrX4EHCTqiC6I0rcdHRUSy8vURHAOcBQwAbsw5vh749zz39ifoWM6qAg7MvcDM+gPjCPoamk0EZnY+QV8Fe+yxR563zU99BCIiW+TrI7gPuM/MvuPuD7byta2JY40b5n9D8ChqvVlTlzfEMRWYClBZWbnDjfs713/G1Pqfw/rhULHLjr6ciEiHFnUcwYNmdgKwH1Cec/yqFm6rAnbP2R8AfNzomkpgRpgEegPHm1mduz8cJa7t9fcr72EEb8IL18GJN+a/QUSkE4u6ZvEUgj6BI4A7gdOAV/PcNgfY28wGAR8BE4Azcy9w90E573Ev8HisSWByX6irYRQE9ZW5dwU/6TK4YnVsbysi0p5FfWroEHf/Z+Bzd/8FcDBb/7W/DXevAy4heBpoKfA/7r7YzC40swt3JOjtNnEBDB3PppKyIMZ0Fxg2HiYuLEo4IiLtQdS5hrJjCL4ys92AtcCgFq4HwN2fAJ5odGxKM9eeFTGW7VfRD8oqSG+updozlNXXQFl39ROISKJFTQSPmVlP4AbgdYJO3zviCipWG1eztP94fvTuSB4+6G3KNqwqdkQiIkWVNxGYWQnwnLt/ATxoZo8D5e7+ZdzBxWLCNOb+dQVL31nM+qPOpaxbWbEjEhEpqrx9BOFiNL/O2a/psEkglEkFxdZYAhGR6J3FT5vZd6ylh/07kGwiqKvXfEMiIlH7CH4I7ATUmVk1wcOX7u7dY4ssRplUkM9qVSMQEYk8oKwi7kAKqVRNQyIiDSI1DZnZc1GOdRTZpiHNQCoikn/20XKCEcW9zWxntp59dLeYY4tN17IUABtrtFyliEi+pqELgEsJvvRfY0siWEew1kCH1LNLKQBffl1b5EhERIov3+yjNwE3mdm/uvstBYopdj26ZgD44qtNRY5ERKT4onYW32JmhwADc+9x99/HFFesenYJEsGXXysRiIhEnX30D8BgYD6QbVh3oEMmgq6lKTIp4wslAhGRyOMIKoEh7t4pRmCZGT26lKppSESE6COLFwH94gyk0Hp0SauzWESE6DWC3sASM3sVqMkedPeTY4mqAHp2LVUfgYgI0RPBpDiDKIaeXTKsXFed/0IRkU4u6lNDL5jZnsDe7v6smXUFUvGGFq8eXTK8uXJ9scMQESm6qFNM/AswE/hdeKg/8HBMMRVEj64ZNQ2JiBC9s/hi4G8JRhTj7m8DfeMKqhB6dillQ02dJp4TkcSLmghq3L3hERszSxOMI+iweoaji9epViAiCRc1EbxgZv8OdDGzo4EHgMfiCyt+PcLRxRpUJiJJFzURXAasARYSTET3BHBFXEEVQna+IfUTiEjSRX18tAtwt7vfAWBmqfDYV3EFFreG+YY0ulhEEi5qjeA5gi/+rC7As20fTuH07BpMRf2FRheLSMJFTQTl7r4huxNud40npMJo6CNQjUBEEi5qIthoZqOyO2Y2Gvg6npAKo3t5mj58ztGvnA3rVxU7HBGRoomaCC4FHjCzF83sReCPwCWxRVUA6VQJPyp/hN3WzYcXrit2OCIiRWNRZ5Y2swywD8FylW+6e1HaVCorK33u3Lk79iKT+0JdzbbH02Vwxeode20RkXbIzF5z98qmzkWtEQCMAfYHRgJnmNk/t0VwRTFxAQwdTzVlwX66CwwbDxMXFjcuEZEiSOQKZVT0g7IKSqmlllJK62ugrDtU7FLsyERECi7WFcrM7FjgJoKZSu9092sbnf8n4Kfh7gbgInd/ozXvsd02rubFHifzh01HcOeQhbBBHcYikkxRE0F2hbJPor5wOOjsNuBooAqYY2aPuvuSnMveA/7e3T83s+OAqcCBUd9jh0yYxtOzFjJv0Uo48eyCvKWISHsU5wplBwDL3f1dADObAYwFGhKBu7+Uc/3LwICI8bSJHl2CqajdHTMr5FuLiLQbca5Q1h/4MGe/ipb/2j8X+FNTJ8zsfOB8gD322GM7Qmlaz64Z6jY7G2vr6VYW9aMQEelcIj015O4vAG8CFeHP0vBYS5r6E7vJPgYzO4IgEfy0qfPuPtXdK929sk+fPlFCjqRnl3Caia80zYSIJFfUFcpOB14FxgOnA6+Y2Wl5bqsCds/ZHwB83MRr7w/cCYx197VR4mkr3TXNhIhI5KahnwFj3H01gJn1IZh0bmYL98wB9jazQcBHwATgzNwLzGwP4CHge+7+Vitj32FanEZEJHoiKMkmgdBa8tQm3L3OzC4BniJ4fPRud19sZheG56cAPwd6AbeHnbV1zY18i0M2EWhxGhFJsqiJ4EkzewqYHu7/I8107OZy9ycIFrHJPTYlZ/s84LyIMbS5LX0ESgQiklyREoG7/8TMTgX+jqATeKq7z4o1sgLYslylOotFJLlaTARm9k1gF3f/i7s/RNCej5kdZmaD3f2dQgQZl/JMCaXpEi1XKSKJlu+pod8A65s4/lV4rkMzM3p2yWi5ShFJtHyJYKC7L2h80N3nAgNjiajAenbNqI9ARBItXyIob+FclxbOdRh7lq7nBx9O1CplIpJY+RLBHDP7l8YHzexc4LV4QiqsM6tnsG/tIq1SJiKJ1eIKZWa2CzALqGXLF38lUAqMc/eVsUfYSJusUAZapUxEEmW7Vyhz91XufgjwC2BF+PMLdz+4GEmgTYWrlG0yrVImIskWdRzBbGB2zLEUVrhKWdprqfYMZfU1mFYpE5EEas2axZ3PxtUs230842qvonr/72uVMhFJpGRPwj9hGm+98TFL357HR3/7Pb7Zt6LYEYmIFFyyawRAz3Caic81lkBEEirxiWCPb3QF4O1VG4ociYhIcSQ+EezZqys9umRYUPVFsUMRESmKxCcCM+PQXes4c8lFGl0sIomU+EQAcG79TIbWLaZu9rXFDkVEpOCS/dRQOLp4JIBByet3w+t3a3SxiCRKsmsE4ehiTwfz59WVlGt0sYgkTrITQTi62OprqCFDyeYa0OhiEUmYZCcCgI2rYfTZXNf/Nh5NH6vRxSKSOMnuIwCYMA2AXrOXc+k73Thi7LfpUeSQREQKSTWC0P4DetCHz+Ge4/UYqYgkihJBaP/+PflBehbdV8/RIjUikihqGgKY3JcedTV8L/tpzL0r+NFjpCKSAKoRQMNjpLVapEZEEkiJABoeI814LdWexuu+hpK0HiMVkURQIsjauJo1+/wTz9aPDvbff6m48YiIFIj6CLKWP0PfuhpOzH4iX7wPk3pAqhSuXFPU0ERE4qQaQVbYT7A5XQ5AffajGfqdIgYlIhI/JYKssJ+gpK4agBSbg+NvTA9qBpP7FjE4EZH4KBHk2rga9j+D+sH/QF340ThAuhzO+B+44yi48x804ExEOpVYE4GZHWtmy8xsuZld1sR5M7Obw/MLzGxUnPHkNWEanDqF1LuzSYc1AgOoq4Y/jIWP5kLVHHjmP4Kk8LvDgsTwycItSSK7vb3nsolm/Uq457jkbEP7iKO9betz0efSXPnbUGydxWaWAm4DjgaqgDlm9qi7L8m57Dhg7/DnQOC34e/iGnwELH+OsD6wrQXTITxrAA+dB2veDM7lbm/vuezI5g9eTs72iTfCC9cXP472tq3PRZ9Lc+VvQ+bezJfdjr6w2cHAJHc/Jty/HMDdf5lzze+A5919eri/DDjc3T9p7nUrKyt97ty5scS8lYcugAUzGr7sG770RUTag1bOfGBmr7l7ZVPn4mwa6g98mLNfFR5r7TWY2flmNtfM5q5ZU6BHOWs3QO99sfDrP5sMsj+5vJnt7T23mS1PLdVT0um3ayjl0/Qu1Fppu4mpPWzrc9Hn0tR2HDMfxJkImvoDuvF3YZRrcPep7l7p7pV9+vRpk+DymjANen8TKs+BgYdBaTes267YzgO3Cdqa2d7ecyWETy2ly0mxudNvl1kdvXvuTCl17Sam9rCtz0WfS1Pb1Ne0+QJacSaCKmD3nP0BwMfbcU3xTJgWtMWd9Rj8+0fw4zdhl6HQc0+o2BV2HhhUzywF3XYNjlkKUmXbd67brjDklCDj99wTzns2+N3Zt0efDdVfBL/bS0ztYVufiz6X5srfxgtoxdlHkAbeAo4CPgLmAGe6++Kca04ALgGOJ+gkvtndD2jpdQvWRyAi0om01EcQ21ND7l5nZpcATwEp4G53X2xmF4bnpwBPECSB5cBXwNlxxSMiIk2Lda4hd3+C4Ms+99iUnG0HLo4zBhERaZlGFouIJJwSgYhIwikRiIgknBKBiEjCxfb4aFzMbA3w/nbe3hv4tA3D6QhU5mRQmZNhR8q8p7s3OSK3wyWCHWFmc5t7jrazUpmTQWVOhrjKrKYhEZGEUyIQEUm4pCWCqcUOoAhU5mRQmZMhljInqo9ARES2lbQagYiINKJEICKScIlJBGZ2rJktM7PlZnZZseOJg5ntbmazzWypmS02s4nh8W+Y2TNm9nb4e+dix9qWzCxlZvPM7PFwv7OXt6eZzTSzN8P/1gcnoMz/Fv6bXmRm082svLOV2czuNrPVZrYo51izZTSzy8Pvs2VmdsyOvHciEoGZpYDbgOOAIcAZZjakuFHFog74kbt/CzgIuDgs52XAc+6+N/BcuN+ZTASW5ux39vLeBDzp7vsCwwnK3mnLbGb9gR8Ale4+lGBa+wl0vjLfCxzb6FiTZQz/v54A7Bfec3v4PbddEpEIgAOA5e7+rrvXAjOAsUWOqc25+yfu/nq4vZ7gC6I/QVnvCy+7DzilKAHGwMwGACcAd+Yc7szl7Q4cBtwF4O617v4FnbjMoTTQJVzwqivBSoadqszu/mfgs0aHmyvjWGCGu9e4+3sEa7q0uKhXS5KSCPoDH+bsV4XHOi0zGwiMBF4BdnH3TyBIFkDfIobW1n4D/F9gc86xzlzevYA1wD1hc9idZrYTnbjM7v4R8CvgA+AT4Et3f5pOXOYczZWxTb/TkpIIGq8bD9Bpn5s1s27Ag8Cl7r6u2PHExcxOBFa7+2vFjqWA0sAo4LfuPhLYSMdvEmlR2C4+FhgE7AbsZGbfLW5URdem32lJSQRVwO45+wMIqpadjpllCJLANHd/KDy8ysx2Dc/vCqwuVnxt7G+Bk81sBUFz35Fmdj+dt7wQ/FuucvdXwv2ZBImhM5f5H4D33H2Nu28CHgIOoXOXOau5Mrbpd1pSEsEcYG8zG2RmpQSdLI8WOaY2Z2ZG0Ha81N1vzDn1KPD9cPv7wCOFji0O7n65uw9w94EE/03/192/SyctL4C7rwQ+NLN9wkNHAUvoxGUmaBI6yMy6hv/GjyLo/+rMZc5qroyPAhPMrMzMBgF7A69u97u4eyJ+gOOBt4B3gJ8VO56Yyvh3BNXDBcD88Od4oBfBEwdvh7+/UexYYyj74cDj4XanLi8wApgb/nd+GNg5AWX+BfAmsAj4A1DW2coMTCfoA9lE8Bf/uS2VEfhZ+H22DDhuR95bU0yIiCRcUpqGRESkGUoEIiIJp0QgIpJwSgQiIgmnRCAiknBKBCIiCadEIIllZr3MbH74s9LMPgq3N5jZ7TG8371m9p6ZXdjCNYea2ZLcqYhF4qZxBCKAmU0CNrj7r2J8j3sJBr3NzHPdwPC6oXHFIpJLNQKRRszs8JxFbiaZ2X1m9rSZrTCzU83sejNbaGZPhnM7YWajzewFM3vNzJ7Kzg+T533GhwutvGFmf467XCLNUSIQyW8wwZoHY4H7gdnuPgz4GjghTAa3AKe5+2jgbuDqCK/7c+AYdx8OnBxL5CIRpIsdgEgH8Cd332RmCwlWx3oyPL4QGAjsAwwFngnmRCNFMGdMPn8B7jWz/yGYUVOkKJQIRPKrAXD3zWa2ybd0rG0m+H/IgMXufnBrXtTdLzSzAwlqG/PNbIS7r23LwEWiUNOQyI5bBvQxs4MhWBPCzPbLd5OZDXb3V9z958CnbD2/vEjBqEYgsoPcvdbMTgNuNrMeBP9f/QZYnOfWG8xsb4IaxXPAG7EGKtIMPT4qUiB6fFTaKzUNiRTOl8B/5htQBjxG0FQkUhCqEYiIJJxqBCIiCadEICKScEoEIiIJp0QgIpJw/x8pUcb7BkhowwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.plot(campsol[:,0],campsol[:,1],label='CAMP')\n",
    "plt.plot(solTchem[:,sp,t_indx],solTchem[:,sp,A_indx],\"*\",label='TChem')\n",
    "plt.xlabel('Time [s]')\n",
    "plt.ylabel('Concentration of A [mol/m3]')\n",
    "plt.legend(loc='best')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7fa870573c10>"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAxPElEQVR4nO3deXxU9dX48c+ZmUwCEkA2aUGbgLuyKMGqrS1qbcW6QdEHq30qal1++itqV9tqcWlr3W2tDw+tSlst+ChuWOtapP09boBYQBSLiBgUAsgaksx2fn/cO8kkTDKTZG5muef9es2LzJ17J+cqzJnvdr6iqhhjjPGvQL4DMMYYk1+WCIwxxucsERhjjM9ZIjDGGJ+zRGCMMT4XyncAnTVo0CCtqqrKdxjGGFNUlixZsllVB6d7regSQVVVFYsXL853GMYYU1RE5MP2XrOuIWOM8TlLBMYY43OWCIwxxueKbozAGGOSotEotbW1NDY25juUglFRUcHw4cMpKyvL+hpLBMaYolVbW0tlZSVVVVWISL7DyTtVZcuWLdTW1lJdXZ31ddY11NbODfDARPhkufPnzo35jsgY047GxkYGDhxoScAlIgwcOLDTLSRLBG0tvAXWvQaPXQQfvgr//SVLBsYUMEsCrXXlv4d1DSXdNARiTS3PN73r/LlrA9x+IATDcO2m/MRmjDEeshZB0vRlxPceQbu7M8QjTrIwxpgUGzZsYOrUqYwcOZJDDz2UU045hffeew+AO++8k4qKCrZv3958/ssvv4yIcN999zUfW7p0KSLCbbfdBsD5559PdXU1Y8eO5cgjj+TVV1/19B4sEYDzAX/7QQS3rqHDRlWsCW5Mu0LbGONDqsqkSZOYMGEC77//PitXruSXv/wlGzc63clz5sxh/PjxPP74462uGzVqFA8//HDz87lz5zJmzJhW59x666289dZb3HzzzVxyySWe3oclAoDpy+Dws4i7/zmUAEgIheYWQnNL4fBv5CFAY0whWrBgAWVlZVx66aXNx8aOHctxxx3H+++/z65du7jpppuYM2dOq+v2228/Ghsb2bhxI6rKs88+y8SJE9P+ji996UusXr3a0/uwMQKAyqFoeR8EJSphyohBv2FIqBe62RkraG4p/GuO8wiVw8/q8hayMaa16+e/zcqPd+T0PQ/9bF9+ftph7b6+YsUKxo0bl/a1OXPmcM4553DcccexatUq6urqGDKkpXt5ypQpPPLIIxxxxBEceeSRlJeXp32f+fPnM2rUqO7dSAbWIgDYuYH4iieYFzuOBcf9BcZNg6GjYND+yOhzaKo6kVhzawEYMAKmL89ryMaYwjZ37lymTp1KIBBg8uTJPPLII61eP/vss3nkkUeaE0ZbP/jBDxg7diyzZs1qNZ7gBWsRACy8hWDTNhoJc9y446DfSa1eLr9+AJAA3JbBp2tsJpExBaajb+5eOeyww3j00Uf3OL5s2TL+/e9/c9JJzmdJJBJhxIgRXH755c3nDB06lLKyMl544QXuvvtuXnnllVbvceuttzJlyhRvb8Dl7xbBTUNgRj9YfB+C8q3Qiwy9c589ZweNPB4GjCAqYcDGC4wxjhNOOIGmpiZ+//vfNx9btGgR06dPZ8aMGaxdu5a1a9fy8ccfs379ej78sHUl6BtuuIFf//rXBIPBng69FX8nAneQWEO9AIhIOYw6a89un/PmQfUEQhoB2owXzOhn00qN8SkR4fHHH+eFF15g5MiRHHbYYcyYMYOXX36ZSZMmtTp30qRJzJ07t9WxY489ljPPPLMHI05PVNudOV+QampqNKcb08y/Cl3yAE0aolxiSM0FcOode54391wIV7Jr6wYq1i0kJE5XEQNGwLRnoXKf3MVkjMnKO++8wyGHHJLvMApOuv8uIrJEVWvSnW9jBPV1LBk8iRs++TyPHfUeoV3tlJOY+hAAfa4fAMkkAC3jBTaLyBhTpCwRTH2In9y5kCFVFYROvzTz+SOPJ7FlDbGtHxEmigKSbBUYY0wR8vcYAbBheyPvbdzFcQcMyu6C8+YR2PYhYaJAm1lENlZgjClCvk8E723cCcDYfftnf5E7iyjWPItIoPcguOjvHkRojDHe8n0i2NbgfLMfsFc4+4vcWURBYjRRBii6ezMsvt+TGI0xxku+TwSRret5OHwDe+vWzl1YX4cglBNFcLuIFt9n00mNMUXH94ngoFUzGS+r6L/ozs5dOPUhuHolHH4W0YBTI8TKTxjjL1u2bGHs2LGMHTuWoUOHMmzYsObnt9xyCwcffDCHH344Y8aM4U9/+hMAVVVVbN68Oc+Rt+bfWUPuRjSjAAQCS+6HJfd3bhpo5VB4+zHKNI77NlZ+wphCt3MDPDoNpszu9vqfgQMH8tZbbwEwY8YM+vTpw/e//31mzpzJ448/zhtvvEHfvn3Zvn07TzzxRLdD94p/WwTuquKIuBX/Qr3SryrOxB04TgRTWgVg5SeMKVTJ7WgX/tqzX/HLX/6Se++9l759+wLQr18/vv3tbze//tvf/pYjjzySUaNG8e67ToXj+vp6LrjgAsaPH88RRxzBk08+CcDs2bM588wzOe2006iuruaee+7hjjvu4IgjjuDoo4/m008/7Xa8/k0ElUOhvJKQRpwB33gTlPft/DcEd+A4EHe2ubTyE8YUqJTaYmjCszG9nTt3snPnTkaOHNnuOYMGDeLNN9/ksssua96V7Be/+AUnnHACixYtYsGCBfzgBz+gvr4ecMpd/+Uvf+GNN97gpz/9Kb1792bp0qUcc8wxzV1O3eHfRABQX8fzvU7h2sF3O6Wn21tVnMX7MPocYiO+YuWqjSlUbi8Abm2xLvcCZKCqGTeQnzx5MgDjxo1j7dq1ADz//PPcfPPNjB07lgkTJtDY2Mi6desAOP7446msrGTw4MH069eP0047DXB2Okte3x3+HSMAmPoQd9y5kBF9+8Cp53brfQBC7ZWrtvITxuSf2wtAvAlCFV3vBcigb9++7LXXXqxZs4YRI0akPSe5CU0wGCQWiwFOApk3bx4HHXRQq3Nff/31VpvWBAKB5ueBQKD5+u7wd4sA2LY7Sv/eZbl5s3Tlqq1VYEzhqK9zWv8Xvdi9XoAMrrnmGi6//HJ27HB2TNuxYwezZs3q8Jqvfe1r/Pa3vyVZCHTp0qWexJaOpy0CETkZuBsIAn9Q1ZvbvN4PeBDYz43lNlV9wMuYUqkq2xqi9MtVIjhvHlw/wGYRGVOo3NY7kL7KcI5cdtll7Nq1i/Hjx1NWVkZZWRnf+973Orzm2muv5corr2T06NGoKlVVVTz99NOexZjKszLUIhIE3gNOAmqBRcA5qroy5ZyfAP1U9UciMhhYBQxVdQv/p5HLMtQNkTiHXPcsPzr5YC6b0P7ATqc8+A34dA2J7esJxJuconQAY86BSTNz8zuMMYCVoW5PZ8tQe9k1dBSwWlXXuB/sc4Ez2pyjQKU4Iyt9gE+B7nd4ZWlbg5Nv+vXKUYsAbBaRMaboeJkIhgEfpTyvdY+lugc4BPgYWA5MV9VEm3MQkYtFZLGILN60KXfdK9vdOkM5GyNIcmcRxUfaLCJjTOHzMhGkmz/Vth/qa8BbwGeBscA9ItJ3j4tUZ6lqjarWDB48OGcBbtvtJoJctgjA6YecPJPgmgWE0s0iujF392CM3xXbLote68p/Dy8TQS2wb8rz4Tjf/FNNAx5Tx2rgA+BgD2NqJZkIcjZY3NYe5apdturYmJyoqKhgy5YtlgxcqsqWLVuoqKjo1HVezhpaBBwgItXAemAq8M0256wDTgT+KSL7AAcBazyMqZXt7hhB/96dKEHdGefNg/lXEfzUKU/darzgX3NsfYEx3TR8+HBqa2vJZZdxsauoqGD48OGdusazRKCqMRG5AngOZ/ro/ar6tohc6r4+E7gRmC0iy3E+J3+kqj1Wlq+5RZDrrqFU9XXI6HNo2L6RsrUvE5KEbW9pTI6UlZVRXV2d7zCKXoeJQEQmZ/Eejar6TLoX3OPPtDk2M+Xnj4GvZvE7PLGtIUooIOwVDnr3S9x5y71SNr23VcfGmEKSqUXwe+BJ0g/8Jn2JNh/2xWJ7g7OqOFNdkJwYeTz66Rqin9YSJmKtAmNMwciUCP6mqhd0dIKIPJjDeHqU7tjAffpz2Dk65/VG9nDePOT6AYSxVcfGmMLS4awhVT0v0xtkc06hOqFuNqPiKz2tS95Kcu+CgO1dYIwpHBlLTIjIUABV3eCWgTgOWKWqb/dAfHvISYkJd3eyPfREf/38q5yd0NKx8QJjjEe6XGJCRC4BXgVeE5HLgKeBU4HHROTCnEfaU9y65I2400Y9qkuelrvqOGGrjo0xBSLTGMEVwGFAL+BDYH+3ZbA3sAC4z+P4vOHWJQ9rlGggTJlHdcnTcmcRBa4fQCDdqmMbLzDG9LBMiSCqqruB3SLyvqpuAFDVrSJS1Ev5ErvqeDB+IsEjLuDc0Eue1SVv18jj4dM1xLfVEkxEWqqU2niBMaaHZSoxkRCR5GqrrycPikhFFtcWtK2n3c91sQuIDTnMqUueWqe8J7hVSoMJZ3WzVSk1xuRLpg/zybhd2Kpam3J8INDxLgsFbltDD6wqzsT2OjbGFIAOu4ZUdV3qc7cyaAhoAN70MC7PeV5wLhu217ExpgBk1b0jIpeIyEZgGbDEfeRmm7A82dHgUQnqrnDXF8QDziymBIGem8VkjPG9bPv5vw8cpqpVqlrtPkZ4GZjXtnldebQzkuMFGnM2vlflk50xePR82NnDg9jGGN/JNhG8D+z2MpCe5tmmNF1VXwfjphG/4EX+Wj6R+Nr/RT98tedWPRtjfCvbMtTXAK+IyOtA85JcVf2uJ1H1gGQi6FsoicAdL6i4aQinpa56Xnyf87DxAmOMR7JtEfw38HfgNVrGCJZ4FVRP2N4QpbIiRDDQA5VHO8Nd9RwPOjsMRaQcDj4Vho62biJjjCeybRHEVPVqTyPpYckS1AXHXfUcTESISphQIsL22nfoV/+B00106h35jtAYU2KybREsEJGLReQzIjIg+fA0Mo9t2x2hf68CGChOxx0vCAUgIEq/Xe+DJpwuIltsZozJsWxbBMm9hq9JOaZA0c4c2laoLQJoHi+QL/+Qxr9eA+8+TQURNNQLOeRU+Oov8hygMaaUZKo++hmAlCmj1aUyfXT77mjhDBS3p3IoFXv1p5wojVqGxhpJSNCmlRpjcipT19D9IvKaiNwsIhNExLPN7nvatoZo4Uwd7Uh9HVJzAQuOm8ODsRPZ/u5CWPeaTSs1xuRMphITE90CcxOAScBtIrIOeBZ4tm0JimKhqoU7WNyW20008aYhEGqCiHvcppUaY3Ik42Cxqjaq6rOqOt3d3eZ7OAnkHhF5w/MIPbCrKUY8oYU7WJzO9GXED5tCE842l4lghZWhMMbkRKdLSavqB6p6r6qeDnzRg5g8VxAF5zqrcijBir6EidJEGcSbaIiLjRcYY7ot02DxThHZ4T52pjzfKSI7VDXS0fWFanshlKDuivo6pGYaH02ez5zEV9jx7kLUxguMMd2UaYygsqcC6UnNLYJiSwTueMH+Nw1h/0BTsnK1jRcYY7ol664hERkjIle4j9FeBuW1hmgcgL3CRToJyi1DEQ044wXRQIWVoTDGdFm2+xFMBx4ChriPh0Tk/3oZmJeiceerdFmowOoMZcstQxHSKBEJE4w3sfPjd2H9EusmMsZ0WrYtgguBz6vqdap6HXA08B3vwvJWMhGEAkW87XJ9HTJuGmVuGYrKHautDIUxpkuy/SQUIJ7yPE7KfuvFJhpXAMLBIk4EUx+CU+9ArlxO48GTacTd3SzUy6aVGmM6JdtPwgeA10VkhojMwClHfZ9nUXkslmwRBIs2l7VoU4aCWCMRDdi0UmNM1rJKBKp6BzAN+BTYCkxT1bs8jMtTzWMExdwiSOWWoVg+cR4Pxb/CtndsWqkxJnudmTbzARBzrxEROVJV3/QmLG+FdtfxcPgGwg2joXLffIfTfe600vE3DWF8sKmlE8+mlRpjspDtrKEbgWXAb4Db3cdtHsblqdHvz2K8rKL3q7fnO5TcsmmlxpguyLZFcDYwslhXEje7aQjEmjgMQCCw9AFY+kDpfGNOmVYadaeV2u5mxphMsu0kXwH07+ybi8jJIrJKRFaLyI/bOWeCiLwlIm+LyMLO/o5OafONWUtxho07rdR2NzPGZCvbFsGvgKUisgJoSh50C8+lJSJB4HfASUAtsEhEnlLVlSnn9AfuBU5W1XUi4u2nVMp+wI1aRkW8Ccr7QuU+nv7aHpWyu1nkbz8hsfJpKmgiEawgcOhptruZMWYP2SaCPwK/BpbTUuEmk6OA1aq6BkBE5gJnACtTzvkm8FhyXwNV9b5/pr6OpUMmc8MnR/Hk+PdgV4n2nVcOJdyrH+pWKy2LN7EzolQ+ej5MmV1ayc8Y0y3ZJoLNqvqbTr73MOCjlOe1wOfbnHMgUCYiLwOVwN2q+qe2byQiFwMXA+y3336dDKONqQ/x9Py3WbOxFk69tHvvVejcaqWbRv4H//vw7Xxx1T/ow2bExguMMSmyTQRLRORXwFO07hrqaPpoutVamub3jwNOBHoBr4rIa6r6XquLVGcBswBqamravkenxeJKWahE1hB0xO0mGn7TEP6j5X+bTSs1xrSSbSI4wv3z6JRjCpzQwTW1QOok/eHAx2nO2ayq9UC9iPwDGAO8h4ei8QShQAmsKs7W9GXw3M+IvzOfYLyRRsoJHfgVQrvrnGml1k1kjK9llQhU9fguvPci4AARqQbWA1NxxgRSPYmz5WUICON0Hd3Zhd/VKdG4ls6q4mykDJLHA+WE4xHWr1nB8PhH1k1kjMm4Q9mpmd6gvXNUNQZcATwHvAP8j6q+LSKXisil7jnvAM/iLFZ7A/iDqq7o3C10XjSeoKwU6gx1Rn0djJtGUJSAKPvGPkRsWqkxhswtgltFZD0dVxr9JfB0uhdU9RngmTbHZrZ5fitwa+ZQcyeWSBDyU4sAmscL+PIP4bmfEVs5n1CikYiUU3bQSciujdZNZIxPZUoEG4FM/Qb/zlEsPSYS81nXUKrm1ccRYoEwoXiEjR+sYJ/IOusmMsanMu1ZPKGH4uhRsYQPu4ZSJbuJlv4ZEWVo01rnuM0mMsaXfPm12Bkj8OWtO1I2tdHDpxCRZJG6ctSK1BnjO778NIzG1V/TR9tTORQp70sZySJ1ETZ/+DZqex8b4ys+TQQJwn5YUJaNNkXqBjd8YLOJjPGZjJ+GIvI5ERnk/ny0iHxfRCZ5H5p3YtYiaGHdRMb4XqZ1BNcCfwdeE5GbgLuAQcB3ReQuz6PziO/HCNJJ0020xbqJjPGFTNNHzwEOAXoD64ChqrrbXQn8lsexecYSQTuS3URL/4yoMqjhA+e4zSYypqRl+jRsVNWIqm4D3lfV3dC8arhodytzSkxY19AerJvIGF/K1CLoLyKTcVYW93V/xn3ez9PIPBSL+3BlcWek6Saq++BthkQ+tEVnxpSgTIlgIXCa+/M/Un5OPi9KEb8VneuKNt1E+zRZN5ExpSrTyuJpPRVIT/L9yuJspGx5qc/9lNjK+ZQlmohIOaEDTyJQb7WJjCkVvvxaHI3ZYHHWkt1EGnVqEyUifPz+cptNZEwJ8eWnYTShhKxFkD23NlFInEVnw62EtTElJdM6grPcP6t7JpyeEY0nCFuLIHvubCKuXA6Hn0UsWAFAI+U07n+KzSYypshl+jS8xv1znteB9JR4QlGFUMASQaclS1gndzpTt5uodrF1ExlTxDLNGtoiIguAahF5qu2Lqnq6N2F5JxpPAFAWsq6hLkkpYY0oI/Qj57jNJjKmaGVKBF8HjgT+DNzufTjea04E1iLomjY7nSXemU8g3kiDhtm93/EM1E9tNpExRSbT9NEITp2hY1V1k4hUOod1V8+El3vRuALY9NHucruJAokIGiynPB7how/fZkBgvS06M6bIZPu1eB8RWQqsAFaKyBIROdzDuDwTc1sEtrI4B9xuIkEJoBwYqEXQltlEM/rbILIxRSDbT8NZwNWq+jlV3Q/4nnus6ETcRGCzhnKgzWwiDfUCoEHDbA4PR8EGkY0pAtl+Gu6lqguST1T1ZWAvTyLyWMztGrJ1BDnkdhNJvAkFekmEQZE2rQNba2BMwco2EawRkWtFpMp9/Az4wMvAvJLY8QkPh2+gT3RLvkMpLcluom89AQNGkHD/ajVSToOtNTCmoGWbCC4ABgOPuY9BQFHWIeq/6C7GyyoOXnVvvkMpLcluopHHQ/UEAkLLWoPVttbAmEImqprvGDqlpqZGFy9e3PkLbxoCsaY9j9u899ybey702QeW/hni6batEPjeKptiakwPEpElqlqT7jX/jJhOXwaHn0XcLY8QD1bAqLNg+vI8B1aCOhhE3lqxrw0iG1Ng/JMIkvPe4000ahmBeBOU97VvpV5KM4i8d+NHNohsTIHxTyIAqK/jkwO+yaTIDdQddC7sssFLz6UMImubQeTdIyfaILIxBSBTiQkARGQw8B2gKvUaVb3Am7A8MvUhVr1bxzvLF/HJF7/JPvv2z3dEpS9ZkgKQ6gnI1tnEpZxwPMKa95czEluJbEy+ZZUIgCeBfwIvAnHvwvFeckFZKGDrCHpcm4J1+1PrHE8WrAuGYXgNTJltXXbG9KBsE0FvVf2Rp5H0kOSCsnDIX71iBaFNwTp992kk1kCDhnm735cY+7lBhFb8jzOQbC0EY3pMtp+GT4vIKZ5G0kOi1iLIv9RB5FA5FRKhZseLhJbPhdSdz6xWkTE9IttEMB0nGTSKyE73scPLwLzSXIbaag3lV3IQ+aKXkNHn0NhrKA2EAXdq74ARznk2zdQYz2XVNaSqlV4H0lNaylBbIsirlEFkJs+kYv5V6JIHUCAYb4RP1ziv2YY3xngu609DETldRG5zH6dmec3JIrJKRFaLyI87OG+8iMRFZEq28XRVLJEsQ21dQwWlvg6puYDINx9jU3g4cXX+aiZCveDgU22aqTEeyioRiMjNON1DK93HdPdYR9cEgd8BE4FDgXNE5NB2zvs18FznQu+aSMy6hgqSuxq5/MATGTzqqwREadQyiDWy6+N3Yf0S6yYyxiPZfhqeApykqver6v3Aye6xjhwFrFbVNe5OZ3OBM9Kc93+BeUCPtPtjCduhrOC5rYNwUAig9Nmx2gaRjfFQZ74W90/5uV8W5w8DPkp5XuseayYiw4BJwMyO3khELhaRxSKyeNOmTdlF246otQgKn9s6CFy1nNhhU4hIOeCsRo70q3bOsdaBMTmT7TqCXwFLRWQBIMCXgGsyXJPuK3fbUqd3AT9S1bhI+9/QVXUW7o5oNTU13SqXGnVbBDZ9tAhUDiVU0ReIokAFTbDd3QYjOYhslUyN6basvhar6hzgaFr2IzhGVedmuKwW2Dfl+XDg4zbn1ABzRWQtMAW4V0TOzCamrorGE5QFhY4SjykgKbWKov2ribt/ZSNSTnxvm2JqTC502CIQkYNV9V0ROdI95NYE4LMi8llVfbODyxcBB4hINbAemAp8M/UEVa1O+V2zgadV9YnO3ULnxOIJQgHrFioaKdNMy0Yej745G1UIaxNsbTPF1FoHxnRJpq6hq4GLgdvTvKbACe1dqKoxEbkCZzZQELhfVd8WkUvd1zscF/BKNK42UFys6uuQcdPgkNNofOJKynauI0iCaKCcYP9hBD79wMpTGNMFHSYCVb3Y/XGiqjamviYiFZneXFWfAZ5pcyxtAlDV8zO9Xy44XUPWIihKKa2DigNPaG4dlCWa9lyAZq0DY7KW7SfiK1keK3iWCEqE2zqQbz1BfZ+q5rGDaKActfIUxnRKpjGCoThTPnuJyBG0zATqC/T2ODZPxOJqq4pLQUrrYK+DTkCXzEax1oExXZHpq/HXgNtwZvzcgTNWcDvO2MFPvA3NG5F4grC1CEpLfR1S47QOdu31OWsdGNNJmcYI/gj8UUS+oarzeigmT1mLoASltA76HHyitQ6M6aRs1xHME5Gvi8gPReS65MPr4LxgYwQlLqV1UN+ndesgYesOjEkr2z2LZ+KMCRwP/AFn8dcbHsblmWhCCVkiKF2txg7atA5s3YExaWX7iXisqv4nsFVVrweOofWq4aIRjSUIW9eQP6S0Dhoqq1qtSo71t5pFxiRlW2souYZgt4h8FtgCVHdwfsGKJWxlsW+ktA56paw7CGsTbLOaRcYkZfuJOF9E+gO3Am8Ca4E5HsXkqUhcKbON6/0nZd1BpF9LzaJGymnqW+WcY60D41MZWwQiEgBeUtVtwDwReRqoUNXtXgfnhVg8QZlVHvWflNZBeP+WmkUVNMGOtc4L1jowPpXxq7GqJkipNaSqTcWaBMBmDRlatQ7i/UeQcP8ZNBBmR+/9nFrp1jowPpLtGMHzIvIN4DFV7dZ+APlm6whMausgOHICuK2DXkTotXud84K1DoyPZPvV+GrgEaBJRHaIyE4R2eFhXN7YuYHb669hoG7NdySmUKTsd6ADUloHGqaubFhL62DnBnhgom2RaUpSVi0CVa30OpAesfAWxug7NGz5I3BivqMxhSCldSDVE5CtbutAIvSKrndeaG4dYGWuTUmSbHp6ROQlVT0x07GeUFNTo4sXL+7cRTcNgVjTnsdD5fCzutwEZorf3HOhzz5wyGnw16vRrWsRTbRzsnUZmeIiIktUtSbdax12DYlIhYgMAAaJyN4iMsB9VAGf9SBWb0xfBoefBaFegLOgiFFnwfTleQ7MFJSpDznf9kceD9UTnFK7wXIUiKVMN93dp8o53waUTYnINEZwCbAEONj9M/l4Evidt6HlUOVQKK+EeBONWkZII1De177Nmfa5Ywd85yVk0MEESaA4001771oLqNNdNKMfzOhvYwemqHWYCFT1bndf4e+r6ghVrXYfY1T1nh6KMTfcf9hnxW7krX0mwy77h2s6kGwdDB0Fg/ZHai5EvvUEib1bDyhvtAFlUwKyGiMAEJFjgSpSBphV9U/ehNW+Lo0RuFSV6mue4bsnHsDVJx2Y48iML8y/yp1umqDdScg1F9qAsik4HY0RZFt99M/ASOAtIO4eVqDHE0F3xBJO0rOVxabLktNNOxpQtjUIpshku6CsBji0FBaTAVZryHRdynRTqp3ppgTL0XgTCYIEidOgYXZX7MOAplpk4a/hyz+ER6fBlNmWFExByvYTcQUw1MtAekIk7nxzC1mLwOTCHgPKcRRnDcLApo+Q5IDy7QfBh6/YLCNTsLJtEQwCVorIG0DzhHxVPd2TqDwScxNB2FoEJhdSWweD9oeqL1iXkSlK2SaCGV4G0VOibteQ7Udgci7LLqP68iEMjKy3LiNTULLds3ghzh4EZe7Pi3D2JSgqUbdFUGZF54yXOugyGhSpTd9lZFNPTR5lO2voO8DFwACc2UPDgJkUWcGelkRgLQLjoS53GWG1jExeZNs1dDlwFPA6gKr+W0SGeBaVR5qnj1oiMD2lgy4jJUigeTa2y8YRTB5k+4nYpKqR5BMRCeGsIygqkZg7a8i6hkw+tOkyChBHg+E9ahltb7s5jnUbGY9l2yJYKCI/AXqJyEnA/wHmexeWN1paBJYITB6k6zKqmQaPXkRw87vNtYwq2m6OI+73NRtgNh7JNhH8GLgQWI5TiO4Z4A9eBeUVGyMwBaNNUpCqL8Ahp6F/vRrccQRVEAGSYwo2lmA8km0i6AXcr6q/BxCRoHtst1eBeSHavKDMEoEpIG02x2HrbAhVQKyRLeXD6d1URy8ira+xsQSTQ9l+Ir6E88Gf1At4MffheCu5jiAcsq4hU6CS4wgXvYjUXMjAXgEqJEZMypyxBHX+yUYD5UT6VTvX2PRT003ZJoIKVd2VfOL+3NubkLwTsxaBKXSp5a9PvQOGjkbGTSN0yQIYdDBBcfZFKEs0Ed7+AdiaBJMD2X4i1ovIkcknIjIOaPAmJO/YGIEpOimJQVL2RYj1r27eF6EVSwqmC7L9RLwSeERE/iki/wQeBq7IdJGInCwiq0RktYj8OM3r54rIMvfxioiM6VT0nZTsGrJZQ6YopWylGRp5PAGhefppnOCe51tSMFnKtsTEIpztKi/DmTp6iKou6egad0D5d8BE4FDgHBE5tM1pHwBfVtXRwI3ArM6F3znWIjAlI7kvwnf+3lLGItB6TUIrVgXVdCDbWUMA42nZoewIEcm0Q9lRwGpVXQMgInOBM4CVyRNU9ZWU818Dhncink5L7kdgC8pM0cuwJiERCCOJCHENEBKrgmo6ltVXY3eHstuAL+IkhPE4m9V0ZBjwUcrzWvdYey4E/pZNPF2V3I8gbC0CU0rS7K8cuNhtKUiCeJsZR02Us2Ov/WyvZdPMyx3K0n3tTnu9iByPkwi+2M7rF+MUvWO//fbrRAitNc8askRgSlWahWrBNquXy2mivL716mUFK43tY17uUFYL7JvyfDjwcduTRGQ0zirlM1R1S7o3UtVZqlqjqjWDBw/uZBgtQrvreDh8A+GGui6/hzFFI01LQb71BDpgBCqt/+kLNI8jqA0u+062iSC5Q9lzIvJU8pHhmkXAASJSLSJhYCrQ6hoR2Q94DPiWqr7X2eA7a8yaWYyXVVS8crvXv8qYwpIy40iqJzgf/MFyFEhI6xlHlhT8R7Lp7RGRL6c77m5S09F1pwB3AUGcEhW/EJFL3WtnisgfgG8AH7qXxFS1w7GHmpoaXbx4ccaYW7lpCMSa9jweKoefWevA+Mzcc6HPPuB2GbH53ZbS2BIkoPG0lykgNRda91GREpEl7X2+ZpUI3DfZB2eQGOANVc3LJ2iXEsHODfDcz4iufIqyRBOEesEhp8JXf2F/kY2/WVLwjY4SQbazhs4G3gDOAs4GXheRKbkL0WOVQ6G8kmAiQqOWQbwJyvvaX1xj2owjUHNhy34JGu9695F1JRWVbGcN/RQYn2wFiMhgnKJzj3oVWM7V17Fk8CR+tfHzPDbuPdhlf0GNaSXN2gRqpiGPXoR00FJoTgqps48A1r1meygUiWzHCJar6qiU5wHgX6nHekqXuoZcP39yBU/+62Peuu6rOY7KmBLWQfcREkTa6T7ag3Ul5VW3u4aAZ90ZQ+eLyPnAX/F48ZcXInG1yqPGdFYH3Ufidh8BqNt9FCfQvHitFat9VLCyrTX0A+C/gdHAGGCWqv7Qy8C8EIsnCFt5CWO6rp2kQDIphCoIknBXNHei9pElhbzqMBGIyP4i8gUAVX1MVa9W1auALSIyskcizKFoPGGrio3JlXRJ4aIXof/nkP6fI3hJsiBegkQgWSW1/aRgg875k2mw+C7gJ2mO73ZfOy3H8XgqmlArQW2MF1IHmq9c1vKzW+YitSCeBsMQj5AgSJD2B51ZeDOC2KBzD8iUCKpUdVnbg6q6WESqvAnJO9FYwkpQG9OTMlRJ7Xgm0v0t1zZXTMVJCqfe4bQWLDHkRKZPxYoOXuvVwWsFKZZQSwTG5EuWaxZaDTq315U0ox96xyEtrQXrTuqWTJ+Ki0TkO20PisiFQIcb0xQiZ4zAuoaMybtMM5GSg86txhecBJGc8S6aAE3sOcaw8Jb0CcK0K1PX0JXA4yJyLi0f/DVAGJjkYVyeiMata8iYgtPOQjbmnosAMvWhlvGFUDnEmtjRez8qGjYS1pYaYs1jDEltu5NSxxhQ61ZK0WEiUNWNwLHufgGHu4f/qqp/9zwyD0TjSkWZJQJjClYHg87N4wuLH6Dv6hegIbrHGEOMACjt7sqWceWzTxNE1kXnCkV3Vhafcc//Y++9wsyedlSOozLG9Kh0q51DFRBrdMYZ3K06285M6lDNhc6fSx6AcdNKbpZSLlYWlwRbWWxMiciwhqF5q07iaDDcqnBehyufF9/Xatwh7YK3Ehx36Mzm9UUvFk8QDtlgsTElJZvupGThvFAFwVgjKpBo02qIqQDiy24lXyWCaDxhLQJj/KITg9DJBAFKXMoIaJQ4AUK0JIWMg9FQtAnCV4mgb2wL0z/6Oex8pGD/hxhjPJCh1UDNNGTuuQAE3QQRSlkFnRyMjhNAoVWCAIo+QfhqsPjR689msj5PoOYCp3/RGGPSyTAYTSAMiQhKkAAdzFZqT0cD0x4liJxsVVkobM9iY0yPSk0KbqsBt9XQ3mylQkwQlgjcPYsblj9JL4nYnsXGmO7roQTRvD80tCSILvRodJQI/DFG4O5ZXE6UqIQpsz2LjTHdlcW4Q9uB6eTMpVCbmUvJBJFuDKLdQeoc9mj4IxEA1NfxUOIrxMf8J+eHX7Y9i40x3uhGgnBmLjktCBKtB6kBgiRa92jkiG8SQeLsB7n2J88wvfIAOOn0fIdjjPGbLBMEtJ8gCFWABz0avkkE0YTT1AqHbB2BMaaAdCJBMPUhWPxAzns0fJMIYnFnUDwUsJXFxpgi0F6C8GDqu2++HkfjTovAylAbY0xrvvlUjLotAtuz2BhjWvNRIrAWgTHGpOObT8VYXBnMVk58bVpJlY81xpju8k0iiMQTfDf0OIM+XdJSCMoYY4xPZg3dNIT9Y03sn7xbD1bmGWNMsfJHi2D6MraOPIMGDTvPQ71g1FkwfXl+4zLGmALgj0RQOZRYqA/lRIkHyj1ZmWeMMcXKH4kACDZs5sH4iSybOM+p3me1howxBvDLGAGwumYGB62dhvYebJvSGGNMCk9bBCJysoisEpHVIvLjNK+LiPzGfX2ZiBzpVSwHvnwpR8m7HLDgUq9+hTHGFCXPWgQiEgR+B5wE1AKLROQpVV2ZctpE4AD38Xngv9w/c8fdnaw/gEDl5qUwo5/NGDLGGJeXLYKjgNWqukZVI8Bc4Iw255wB/EkdrwH9ReQzOY0i3RaVHR03xhif8TIRDAM+Snle6x7r7DmIyMUislhEFm/atKlzURz/c8DZ7i31T066sXPvY4wxJcrLRJCuulvbDZKzOQdVnaWqNapaM3jw4M5F8eWrW/2i5l/4he927n2MMaZEeZkIaoF9U54PBz7uwjm5M/Ikz97aGGOKlZfTRxcBB4hINbAemAp8s805TwFXiMhcnEHi7ar6Sc4jmbE9529pjDGlwrNEoKoxEbkCeA4IAver6tsicqn7+kzgGeAUYDWwG5jmVTzGGGPS83RBmao+g/Nhn3psZsrPClzuZQzGGGM65psSE8YYY9KzRGCMMT5nicAYY3zOEoExxvicOOO1xUNENgEfdvHyQcDmHIZTDOye/cHu2R+6c8+fU9W0K3KLLhF0h4gsVtWafMfRk+ye/cHu2R+8umfrGjLGGJ+zRGCMMT7nt0QwK98B5IHdsz/YPfuDJ/fsqzECY4wxe/Jbi8AYY0wblgiMMcbnfJMIRORkEVklIqtF5Mf5jscLIrKviCwQkXdE5G0Rme4eHyAiL4jIv90/9853rLkkIkERWSoiT7vPS/1++4vIoyLyrvv/+hgf3PNV7t/pFSIyR0QqSu2eReR+EakTkRUpx9q9RxG5xv08WyUiX+vO7/ZFIhCRIPA7YCJwKHCOiBya36g8EQO+p6qHAEcDl7v3+WPgJVU9AHjJfV5KpgPvpDwv9fu9G3hWVQ8GxuDce8nes4gMA74L1Kjq4Thl7adSevc8Gzi5zbG09+j+u54KHOZec6/7OdclvkgEwFHAalVdo6oRYC5wRp5jyjlV/URV33R/3onzATEM517/6J72R+DMvAToAREZDnwd+EPK4VK+377Al4D7AFQ1oqrbKOF7doWAXiISAnrj7GRYUvesqv8APm1zuL17PAOYq6pNqvoBzp4uR3X1d/slEQwDPkp5XuseK1kiUgUcAbwO7JPc+c39c0geQ8u1u4AfAomUY6V8vyOATcADbnfYH0RkL0r4nlV1PXAbsA74BGcnw+cp4XtO0d495vQzzS+JQNIcK9l5syLSB5gHXKmqO/Idj1dE5FSgTlWX5DuWHhQCjgT+S1WPAOop/i6RDrn94mcA1cBngb1E5Lz8RpV3Of1M80siqAX2TXk+HKdpWXJEpAwnCTykqo+5hzeKyGfc1z8D1OUrvhz7AnC6iKzF6e47QUQepHTvF5y/y7Wq+rr7/FGcxFDK9/wV4ANV3aSqUeAx4FhK+56T2rvHnH6m+SURLAIOEJFqEQnjDLI8leeYck5EBKfv+B1VvSPlpaeAb7s/fxt4sqdj84KqXqOqw1W1Cuf/6d9V9TxK9H4BVHUD8JGIHOQeOhFYSQnfM06X0NEi0tv9O34izvhXKd9zUnv3+BQwVUTKRaQaOAB4o8u/RVV98QBOAd4D3gd+mu94PLrHL+I0D5cBb7mPU4CBODMO/u3+OSDfsXpw7xOAp92fS/p+gbHAYvf/8xPA3j645+uBd4EVwJ+B8lK7Z2AOzhhIFOcb/4Ud3SPwU/fzbBUwsTu/20pMGGOMz/mla8gYY0w7LBEYY4zPWSIwxhifs0RgjDE+Z4nAGGN8zhKBMcb4nCUC41siMlBE3nIfG0RkvfvzLhG514PfN1tEPhCRSzs45zgRWZlaitgYr9k6AmMAEZkB7FLV2zz8HbNxFr09muG8Kve8w72KxZhU1iIwpg0RmZCyyc0MEfmjiDwvImtFZLKI3CIiy0XkWbe2EyIyTkQWisgSEXkuWR8mw+85y91o5V8i8g+v78uY9lgiMCazkTh7HpwBPAgsUNVRQAPwdTcZ/BaYoqrjgPuBX2TxvtcBX1PVMcDpnkRuTBZC+Q7AmCLwN1WNishynN2xnnWPLweqgIOAw4EXnJpoBHFqxmTyv8BsEfkfnIqaxuSFJQJjMmsCUNWEiES1ZWAtgfNvSIC3VfWYzrypql4qIp/HaW28JSJjVXVLLgM3JhvWNWRM960CBovIMeDsCSEih2W6SERGqurrqnodsJnW9eWN6THWIjCmm1Q1IiJTgN+ISD+cf1d3AW9nuPRWETkAp0XxEvAvTwM1ph02fdSYHmLTR02hsq4hY3rOduDGTAvKgPk4XUXG9AhrERhjjM9Zi8AYY3zOEoExxvicJQJjjPE5SwTGGONz/x+yvqoE7K9HPQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.plot(campsol[:,0],campsol[:,3],label='CAMP')\n",
    "plt.plot(solTchem[:,sp,t_indx],solTchem[:,sp,B_indx],\"*\",label='TChem')\n",
    "plt.xlabel('Time [s]')\n",
    "plt.ylabel('Concentration of B [mol/m3]')\n",
    "plt.legend(loc='best')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7fa840389a60>"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAo6klEQVR4nO3deZhU5Zn38e9d1U13o4CyK0ga0IgSXFtNzGZc4hKNozEJOpmocTTm0qRJTN7EzGSGLBoQdeI4MQ5BNDF5ZcYFo45xiQLJvBoFxLSIG4pLgyyiSAPS6/3+cU411UV112no011d5/e5rr6oc+pU1X2Afu56dnN3REQkuVJ9HYCIiPQtJQIRkYRTIhARSTglAhGRhFMiEBFJuLK+DqC7hg8f7tXV1X0dhohIv7J06dJ33H1Evuf6XSKorq5myZIlfR2GiEi/YmZvdPacmoZERBJOiUBEJOGUCEREEq7f9RHk09zcTH19Pdu3b+/rUIpKZWUlY8eOpby8vK9DEZEiVhKJoL6+nkGDBlFdXY2Z9XU4RcHd2bhxI/X19YwfP76vwxGRIlYSTUPbt29n2LBhSgJZzIxhw4apliRSKhrWwq2nQsO6Hn/rkkgEgJJAHvo7Eeknsgv5zh4vugbe/CssmtnjH18STUMiIv1Cw1q460I45zbAdzzOLeTDx82tbZS98QR23Yd3vMeSW4Kfsgr45/U9ElbJ1AiKwdq1a5k6dSoTJ07k4IMP5rTTTuPll18G4N/+7d+orKzk/fffb79+4cKFmBm33HJL+7lly5ZhZlx77bUAXHDBBYwfP57DDjuMI444gieffLJ3b0pEui/Kt/pF18AbT8B1Hw4Kdm/bUciHj8uX3cpO9fqyKpjyRah9rsfCVSLoIe7OWWedxXHHHcerr77KihUruPrqq1m3LmjPu+OOOzjqqKOYP39+h9dNmTKF//qv/2o/njdvHoceemiHa2bNmsWzzz7LjBkz+PrXvx7/zYhINLtb4OdocaMlLJZbSdEWPnbSOEC6AloboWIwDBrVY7ehpqEesmDBAsrLy7n00kvbzx122GEAvPrqq2zZsoVZs2Zx9dVXc8EFF7RfM27cODZv3sy6desYOXIkDz30EKeddlrez/jUpz7FypUr47wNESkku3knt0knU+Bn5CnsISjwzYw0bbSSIkUbnhpAuq0Jw6GsknRLONCjrBJr2Q4jJsEX5sCSW2FLz3YYl1wi+PH9z7NizeYefc+D9x3Mv54xuctrli9fzpFHHpn3uTvuuINzzz2XT37yk7z00kusX7+ekSNHtj9/zjnncOedd3L44YdzxBFHUFFRkfd97r//fqZMmbLrNyIi0XXVnr9LBX6aFK14egDp1pwCf8Qk7AtzYN7fBy+c+vuOjzOF/+gpcPr1PX6rJZcIitG8efOYP38+qVSKs88+mzvvvJPLLrus/fkvfelLfPnLX+bFF1/k3HPP5Yknnujw+u9973v87Gc/Y8SIER36E0SkB0TowN3W3EpVbqdtlnwFPlEL/OxCflrdjjfNfhxD4Z+t5BJBoW/ucZk8eTJ33XXXTufr6up45ZVXOOmkkwBoampiwoQJHRLB6NGjKS8v59FHH+WGG27YKRHMmjWLc845J94bEEmSTpp3HIc8o3QG5nmL9iadTgp8ohb4MRfyUaizuIccf/zxNDY28utf/7r93OLFi6mtrWX69Om8/vrrvP7666xZs4bVq1fzxhsdV4T9yU9+wsyZM0mn070dukjp6qIz1/N04NqSuTuN0unQaWvh72e6gjRt2IhJpC5+HNvrQ7DXh+Af/wQ1F8Gw/XcU+NPqdjTpTP19795/RCVXI+grZsb8+fOZNm0aM2bMoLKykurqahYuXMivfvWrDteeddZZzJs3j2OOOab93LHHHtvbIYuUpk6+7W9r6ti8k1vgt5AiBcG3fNJAK5auIN3aGFzQWadtEX7D7y5z976OoVtqamo8d2OaF154gYMOOqiPIipu+ruRktVp2/5MfMktO4+/z6ONNEYrpCuwrAKfQs07RfrNvitmttTda/I9pxqBiPQP2QX/oFEdvu1vaWxhjy6+7beSwgi/7Vsa86DwT7U29rv2/DgoEYhI8crTzOPXHxQU5BlLbmHPPC/NfNtPSvPO7lBnsYj0vW506maSQKZVu7PO3BSt2IhJcPFjQUduvs7cIu7A7U2qEYhI3+ikU7fVnVQXzTxNVkFj1Sj23PYWXjaAdIu+7e8uJQIR6T25hX+eGbr5BlBnN/MMaGtmwACDg78GNRf26gzcUqVEICLxKlT4Z+kwhLOzTt1MYZ8p6PVtf7cpEfSAjRs3csIJJwDBUtTpdJoRI0YAcN555zF37lzKyspIp9NcccUVfPWrX6W6upolS5YwfPjwvgxdJB45hb93sTxDZoYu6QGUtTYFJ7tq5lFh3+OSmwhyh6LthmHDhvHss88CMH36dPbcc0+++93vcvPNNzN//nyefvppBg8ezPvvv8+9996726GLFKWs3ylfNLPDUg257fzZ4/fTUYZwqvCPVXITQfbysTH9J7v66qtZsGABgwcPBmDIkCGcf/757c/feOON3H///TQ3N3PnnXcyadIktm7dyje/+U2ee+45WlpamD59OmeeeSa33XYb9957L62trSxfvpwrrriCpqYmbr/9dioqKnjwwQcZOnRoLPch0kEnE7laFswkHRb++SZzZbfz523qUadun0leIvjZSMiMMoBYtn0DaGhooKGhgYkTJ3Z6zfDhw3nmmWe46aabuPbaa5kzZw5XXXUVxx9/PHPnzmXTpk0cffTRnHjiiUCw1PWyZcvYvn07+++/PzNnzmTZsmV8+9vf5re//S3Tpk3rsfhFOuhs2YasVTnzFSYF2/n1bb8oJC8R1NbBw/8MLz4ALR8E274ddDp89qoe/Rh3L7h5/Nlnnw3AkUceyT333APAI488wn333de+VeX27dt58803AfjMZz7DoEGDGDRoEEOGDOGMM84Agl3O6urq8nyCyG6IMMIn36qc2YW/qfDvF5KXCAaNhopBwXZvZZWxbPsGMHjwYPbYYw9ee+01JkyYkPeazAY06XSalpYWIEggd999NwceeGCHa5966qkOG9akUqn241Qq1f56kd3SAyN8VPj3P8mcWbx1PRx5YTDL8MgLe3zbt4wrr7ySyy67jM2bgx3TNm/ezOzZs7t8zcknn8yNN95IZjHAZcuWxRKbSLvcvXazZ/LmaCGFA54eQBltwWifssogCWRm8Wrmbr+TvBoBdPyPGeO3lG984xts2bKFo446ivLycsrLy7niiiu6fM2PfvQjpk2bxiGHHIK7U11dzQMPPBBbjCJRvvmnaYN0BWUa4VOSYl2G2sxOAW4A0sAcd5+R8/wQ4HfAOIKkdK2739rVe2oZ6u7R343k1bAWrp+0Y8GeHNmFf96mHn3L73e6WoY6tqYhM0sDvwROBQ4GzjWzg3MuuwxY4e6HAscB15nZgLhiEkm0rCag9x+6ijZ33rJ9+MCDX7kdzT4VlIW7b5maehIhzqaho4GV7v4agJnNA84EVmRd48AgC4bX7Am8C6jXU6SnZHX+fvCnn1MZDvUcEj69n78NFvwiltGmTt6EijMRjAHeyjquB47JueY/gPuANcAg4Mvu3pb7RmZ2CXAJwLhx4/J+WJThmknT33afkx6SVfhnT/KqynetpWDv8djnrocX7lPhn1BxJoJ8pXJuyXQy8CxwPDAReNTM/uLumzu8yH02MBuCPoLcN62srGTjxo0MGzZMySDk7mzcuJHKysq+DkV6QyfLO+T9Bbc0hEM9aWuGCcfBxPBHEinORFAP7Jd1PJbgm3+2C4EZHnx1XWlmq4BJwNPd+aCxY8dSX1/Phg0bdifeklNZWcnYsWP7OgyJSycLu+X9BpY1zp98nb+SaHEmgsXAAWY2HlgNTAXOy7nmTeAE4C9mNgo4EHitux9UXl7O+PHjdzNckX4mZ1XP3ASgSV4SVZeJwMzOjvAe2939wdyT7t5iZpcDDxMMH53r7s+b2aXh8zcDPwVuM7PnCP4ff9/d3+nuTYgkRs6wz50K/6xVPVX4S1SFagS/Bv5A/vb+jE8BOyUCgDBBPJhz7uasx2uAz0aKVCSpwiYgP+dW1j/wU0a480bbKEbbe1RZU4dVPVX4y64olAj+6O5f6+oCM/tdD8YjIjl7ZTQ/PoOyN57ArjuQzIpY41M72vVTtKrwl93SZSJw968UeoMo14hIAXmWefbrD8K8lfJ814fDPtGwT+kBBTuLzWw0gLuvNbMRwCeBl9z9+biDE0mMPOv9ZNpjPXzcYeSPhn1KDyrUWfx14AfBQ5sJXAA8D/zczK5x952XJxSRaAqs99OcqsAG70vZptehrCL/Hr4iPaBQjeByYDJQBbwB7B/WDPYGFgBKBCLdkTv232FT5X5Ubl9HFU0dFnsrb2sG2qDma1Bzodr/JTaFEkGzu28DtpnZq+6+FsDd3zMzrV8g0l05TUAG7L09WInFgXS42FuHb/2ZQl+Fv8SkUCJoM7Nyd28GPpc5aWaVJHVTG5HuKtAE5KSwoVrvR/pOoURwNuH6QO5en3V+GND1DisiSbZTE5CzvmwMg5s3UGVNtJIiRVsw9l8dv9LHCg0ffTP72MwGh6/5AHgmxrhE+rc8TUCjWla3L/mczrfks0gfibTWUDh66CcECSBTv3Ug/67sIklUqAnIUpiWfJYiFHXRue8Ck7UOkEiOPE1Aa9P7slfLO2ETUJpUZvkHNQFJkYqaCF4FtsUZiEi/lKcJaJ/WNVlNQHmWfxApMlETwZXAE2b2FNCYOenu34olKpFipiYgKTFRE8F/Ao8DzwE7bSUpkihhE9CG8jEMatIoIOn/oiaCFnf/TqyRiBSzPPsAjGzWKCApDVETwYJwA/n76dg09G4sUYkUg9xN4N15w0czinepsiY1AUnJiJoIMltMXpl1TsNHpbRlbQWZ+UWptrXtTxuoCUhKQqHVR/dx97fdXRsCSzIU2Aoy7z4AIv1coRrB3HCl0YXAQ8D/untL7FGJ9JWwI3hd2RiGhMtBZG8FqX0ApBQVWmLi1HCBueOAs4BrzexNgqTwUO4SFCL9Up5awOis5SDybgUpUkIK9hG4+3bCgh/AzMYDpwL/YWaj3f3oeEMUiZcvnAnuvMloRro6giV5onYWt3P3VcBNwE1mNqDnQxLpBVm1gEw/wIdY294poI5gSZJCncUN7FhkbqctVN19cIyxicSm8bEZlLvzetso9km9RxU7agHqCJakKdRHMKi3AhGJXcNa/PpJmDsV4akJqR2FvWoBklSRdxkzs0PN7PLw55A4gxLpMQ1r4dZToWEdmx++CnfntbZRbM+kAkvB0InwD3+AIy9ULUASKep+BLXAxcA94anfm9lsd78xtshEekLWpLBMO2Z2LSA4cZxqAZJoUTuLLwKOcfetAGY2E3gSUCKQ4qNJYSLdEjURGNCaddxKnt8vkaKw6Brc4d0BYxnYuL59UliKVtCkMJGdRE0EtwJPmdn88PjvgFtiiUhkV+SpBQxrqtekMJEIIiUCd7/ezBYCnyD4HbvQ3ZfFGZhIt4S1gHcGjGHPxszSEClSQzUpTKSQ7kwoWwW0hK8xMzvC3Z+JJyyRiH42ElqCldENGNGUtTSEoSYgkQiijhr6KXABwd7FmQlmDhwfT1giETSspXXkFF7aOpDxm/7aeS1ARLoUtUbwJWCiuzfFGYxId2z4n58ybM1S0m37UpFqxssqSLWqI1iku6ImguXAXsD6+EIRiaBhLX7dJAxnRHjqwNTq4IG3aVKYyC6Imgh+Diwzs+V03Kry8129yMxOAW4A0sAcd5+R55rjgF8A5cA77v7piDFJAm1++Cr2xFnVNoqx6U0M8EYoq4KDTofPXgWDRvV1iCL9TtRE8BtgJvAc0BblBWaWBn4JnATUA4vN7D53X5F1zV4EK5me4u5vmtnIbsQuSZG1RlCH2cGZ3qrWRqgYrCQgsouiJoJ33P3fu/neRwMr3f01ADObB5wJrMi65jzgnswGN+6upifZSWal0FVtoxiT2kQFjZodLNKDoiaCpWb2c+A+OjYNdTV8dAzwVtZxPXBMzjUfBsrDOQqDgBvc/be5b2RmlwCXAIwbNy5iyNKvZU0Qy7dSaHDiOHUKi/SAqIng8PDPj2adKzR8NN8SFJ5zXAYcCZwAVAFPmtlf3f3lDi9ynw3MBqipqcl9DylBmV3DMvsFVNKkWoBITKLOLP7MLrx3PbBf1vFYYE2ea94JF7PbamZ/Bg4FXkaSJ8+uYeNVCxCJXZf7EZjZ6YXeoItrFgMHmNn4cEvLqQRNS9n+AHzSzMrMbCBB09ELhcOWkrToGjzcO/gDD3ZBde0XIBK7QjWCWWa2mq5XGr0aeCD3pLu3mNnlwMMEw0fnuvvzZnZp+PzN7v6CmT0E1BGMRprj7st35UakH8tZJmKc9g4W6VWFEsE6oNAKXa909oS7Pwg8mHPu5pzjWcCsAp8hpahhLdx1IVv/4WGW3/kTDmn4X6pMeweL9LZCexYf10txSBItugZ/40n+dPsMtjSmODodLBNhWiZCpFd1Z/VRkd2XZ9+AM1seChoPLQX/+Jj2CxDpZUoE0rvCfQPeqxhL1fZg9zAvq8Kyl4jQfgEivUqJQHpHTofw0MYdu4eZlogQ6VMFE4GZfQjY6u7vmNlHCXYpe9Xd5xd4qUigYS2MmsKWAcMpW/U4lWjfAJFi0mUiMLMfEWxI4+FaQScCC4HPmdmn3X1a3AFKCVh0Db56KW/7GCZaM23pClLaQF6kaBSqEZwLHAQMBN4ERrv7NjMrA56NOTbpz/J0Ch9g9eGT2jdApJgUSgTbw13JmszsVXffBu2TxbRbmXQuz2byO3UKi0hRKJQI9jKzswm+0A0OHxMeD4k1MumfuthMXp3CIsWpUCJYBJwRPv5z1uPMscgOYafw1ooRpF97nEoa1Sks0g8Umll8YW8FIiUg0ynMGCbQpE5hkX5C8whk9+U0B+2POoVF+hMlAtk9YXPQG02DGbn+f9UpLNIPFdqP4Ivhn+N7Jxzpb3zhTNpWL2X7upeosMyiceoUFulPCtUIrgTuBO4Gjog/HOkXcnYSM+DA1OrgOVdzkEh/UygRbDSzBcB4M8vdXQx3/3w8YUlRC+cIrCsbw5BmzREQ6e8KJYLPEdQEbgeuiz8cKWo5ncKjW1bv2ElMzUEi/Vah4aNNwF/N7Fh332Bmg4LTvqV3wpOiEXYKb68cCa8+1mGOgHYSE+nfoo4aGmVmjwBDATOzDcD52l84QcI5AmsYS7U30VamOQIipSJqIpgNfMfdFwCY2XHhuWPjCUuKRk5z0ATeCpuD1CksUiqiJoI9MkkAwN0XmtkeMcUkxSJsDtpUNoyKNxZQRRNtZZWkDjpDncIiJSRqIngt3Jvg9vD4K8CqeEKSohE2B20I9xHwdAWp1iZ1CouUmKiJ4GvAj4F7wuM/A1qHqFTlNAdpHwGR0hYpEbj7e8C3Yo5FikHYHLTBhrLnWws1R0AkAbTWkHQUNge917Yvw1JBc5DmCIiUNiUCCeQ0B304s2SEmoNESp4SgbQ3B61nbwbVL1JzkEjCREoEZjYCuBiozn6Nu38tnrCkV4XNQZva9mW4moNEEidqjeAPwF+APwGt8YUjvUrNQSJC9EQw0N2/H2sk0rvC5qB3UkPZ402NDhJJsi43psnygJmdFmsk0rsyo4PefD7YUEbNQSKJFbVGUAv80MyagObwnLv74HjCktjsNFlMzUEiSRd1QtmguAORXlJbx/t/+D4DVj4Yrh1URUrNQSKJFnn4qJl9HvhUeLjQ3R+IJySJTcNatt9xPn9ZU8FpZNYOUnOQSNJF6iMwsxkEzUMrwp/a8Fyh151iZi+Z2Uoz+0EX1x1lZq1mdk7UwKX7tj56NQPWPMVhvMDmyV/FLn5MzUEigrl74YvM6oDD3L0tPE4Dy9z9kC5ekwZeBk4C6oHFwLnuviLPdY8C24G57n5XV7HU1NT4kiVLCsYsWbL6BTooq4B/Xt/78YhIrzOzpe5ek++5qKOGAPbKejwkwvVHAyvd/bVwy8t5wJl5rvsmcDegEikODWtpGTmFJwZ8jA98QHCurAqmfBFqn+vb2ESkKETtI/g5sMzMFhAMNvkUcGWB14wB3so6rgeOyb7AzMYAZwHHA0d19kZmdglwCcC4ceMihiwALQtnklqzlOFt+1KZaoaySlC/gIhkiTpq6A4zW0hQWBvwfXdfW+Bllu+tco5/Eb5Xq1m+y9s/fzbB1pjU1NQUbsuS9uagzD9w+6xhb1W/gIh00GUiMLNJ7v6imR0RnsrsULKvme3r7s908fJ6YL+s47HAmpxraoB5YRIYDpxmZi3ufm/UG5A8Gtbio6bwwpY9GL/pSaqsKWgO0jBREcmjUI3gOwRNMtflec4JmnQ6sxg4wMzGA6uBqcB5Hd7AfXzmsZndBjygJNADFl0Dq5eSbtuXCjUHiUgBXSYCd78kfHiqu2/Pfs7MKgu8tsXMLgceBtIEI4KeN7NLw+dv3vWwJa+cWcMHqjlIRCKIOnz0GXc/otC53qDho11oWMuGu7/Hnqse0iJyItJBV8NHC/URjCYY/VNlZoezowN4MDCwR6OU3dOwlm3/96s8saaCM7SInIh0Q6E+gpOBCwg6eq/POt8A/DCmmGQXbHv051S+/TRH2gi2HXI+ex77j7DkVjUHiUhBUZuGvuDud/dCPAWpaSiHZg2LSAS73DSU4e53m9nngMlAZdb5n/RMiLJLwmGiz7xXycFbn955mKiISARR9yy+maBP4DPAHOAc4OkY45Iows1l9tSsYRHZDVGXmDjW3Q8xszp3/7GZXQfcE2dg0oWs5qAUGiYqIrsnaiLIzCHYZmb7AhuB8V1cL3GqrdMwURHpMVFXH73fzPYCZgHPAK8Dd8QUk3Ql3Fzmydc3aa9hEekRBWsEZpYCHnP3TcDdZvYAUOnu78cdnOys6fEZDFjzFEcwgi1Tvsrgj1+sYaIislsKJgJ3bwv7BD4WHjcCecYrSqzCfoFwRwHGsh6e+w28ME/DREVkt0RtGnrEzL5gXa0VLfGqrWPFsJO1uYyI9LioncXfAfYAWsxsO8FSE+7ug2OLTHZoWMvG2/6eFeuqmFTWjJepX0BEek7UCWWD4g5EOvfuH3/G3u8s5RMDRtJ2+NcoO+pC9QuISI+JOqHsMXc/odA56WFhv8BQAIPRbetg6S3wt9+pX0BEekyXfQRmVmlmQ4HhZra3mQ0Nf6qBfXslwqQKl49YXPVx9QuISKwK1Qi+DkwjKPSXsmMZ6s3AL+MLSzLLRwzW8hEiErNCO5TdANxgZt909xt7KaZk0/IRItLLonYW32hmxwLV2a9x99/GFFdy1dax5f4fkH75f6hCy0eISPyidhbfDkwEngVaw9MOKBH0pIa1tN15IU++XcEJaPkIEekdUecR1AAHe5RdbGTXLboG3nySD7cNZ82Hz2Psid/QMFERiV3URLAcGA28HWMsyZXTL/Ch1AZY+Xt4/S4NExWR2EVdYmI4sMLMHjaz+zI/cQaWKLV1NA6uJlPfcg0TFZFeFLVGMD3OIBLtpyOhtZEKaB+cay0fwPK74Qtz+jIyEUmISDUCd19EsAdBefh4McG+BLKbfPJZONCW6X1JV8DQiTBRk7ZFpHdEHTV0MXAJMJRg9NAY4GZApdWuCvsFMjP02td1bW2ECcfB6df3UWAikjRR+wguAz5OMKMYd38FGBlXUIlQW0fTkKx+AUvD/ifBoedplJCI9KqofQSN7t6U2Y7AzMoI5hHIrgj7BQbAjn4Bb4VXH4N/fa8vIxORBIpaI1hkZj8EqszsJOBO4P74wipxk89Wv4CIFI2oNYIfABcBzxEsRPcgoCEt3ZU1X8BQv4CIFIeoiaAKmOvuvwYws3R4bltcgZWk2jqabzmFsk2rgi3eLI1NPB72GKF+ARHpM1ETwWPAicCW8LgKeAQ4No6gSlJYGyjPOqV+AREpBlH7CCrdPZMECB8PjCekElVbx6aq/dpHCalfQESKRdQawVYzO8LdnwEwsyOBD+ILq8SEo4T2gh1b+7Q2wnur4FualycifStqIpgG3Glma8LjfYAvxxJRCWo5+CzSz83DCatg6QoYMhaGTujjyEREom9Ms9jMJgEHEnynfdHdmwu9zsxOAW4A0sAcd5+R8/zfA98PD7cA33D3v3Uj/uIW9gtk/pIzlQGNEhKRYhK1RgBwFDt2KDvczLrcoSwcWfRL4CSgHlhsZve5+4qsy1YBn3b398zsVGA2cEw376F41dax7T8/S1XDG8FQUUuDRgmJSJGJc4eyo4GV7v5a+B7zgDOB9kTg7k9kXf9XYGzUwIte2C8wEHZUBTRKSESKUJw7lI0B3so6rqfrb/sXAX/M94SZXUKw6B3jxo3rRgh9aPLZeN0duEPKUL+AiBStOHcoszzn8iYSM/sMQSL4RL7n3X02QbMRNTU1xb3GkWYPi0g/EzURZHYoexpozJx098938Zp6YL+s47HAmtyLzOwQguUqTnX3jRHjKV61dbTccgppzR4WkX4izh3KFgMHmNl4YDUwFTgv+wIzGwfcA/yDu7+8C59RXHJGCYFmD4tI8evODmUvAoPCnxfCc129pgW4HHgYeAH4b3d/3swuNbNLw8v+BRgG3GRmz5rZkl28j+JQW0fDwHGaPSwi/YpF6f81sy8Bs4CFBE3fnwS+5+53xRpdHjU1Nb5kSRHmi3CU0E4spdqAiPQ5M1vq7jX5novaNPRPwFHuvj58wxHAn4BeTwTFyiefBXXzNEpIRPqdqIkglUkCoY1EX7CutGnvYRHp56ImgofM7GHgjvD4y3Qy5j9xNEpIRPq5qGsNfc/MziYY52/AbHefH2tk/UHYL6BRQiLSn3XZvGNm+5vZxwHc/R53/467fxvYaGYTeyXCYqa9h0WkBBSqEfwC+GGe89vC587o4Xj6B80eFpESUqjDt9rd63JPuvsSgpVIk6m2DiadTlvYReykYf+T4NDz1C8gIv1OoRpBZRfPVfVkIP3KoNFse/slqtxpsxQpc9hrnGoCItIvFaoRLDazi3NPmtlFwNJ4QipyPxsJ04cw8P1XMIMUbeBtsPTWvo5MRGSXFKoRTAPmhzuJZQr+GmAAcFaMcRWv2jre/eWJ7P3BW0HfQFkVHHQ6fPaqvo5MRGSXdJkI3H0dcGy4TPRHwtP/4+6Pxx5ZMQqHiw6FHYtst3wAy++GL8zpw8BERHZd1HkEC4AFMcdS9NqXkUCb0ItI6ejOnsXJlbuMROa8houKSAnQekFR1NbRstf49u3V3DRcVERKh2oEhWizGREpcaoRFHLRn2iqGEqrhw1CWkZCREqMEkEBrUtupbzxXVLmeFkFtDUH/QJf0VYMIlIa1DTUmbBJKJ19rqUx2HFM/QIiUkJUI+hMbR3NQ8bv2H+4rAqmfBG+8yJM/X2fhiYi0pNUI8gnrA2UgyaOiUjJU40gn9o6tu35oR21AXUQi0gJU40gV1gbGAg7agOtjfDeKvjWM30YmIhIPFQjyHXRn9hatreGi4pIYigR5Gh4Yg4Dm98L9hgoq9RwUREpeWoaygibhAZBVgfxdg0XFZGSpxpBRm0dW7M7iDVcVEQSQjUCaN9nYA/QcFERSRzVCAAmn40DbRouKiIJlOwaQdgvAEFFwLKHi2qfARFJiGTXCGrr8L0naJ8BEUm05NYIcnYdA+0zICLJlNwaQW0djYOr20cJufoFRCShklkjCEcJVUD7KCHTMhIiklDJSwQNa6FiEL6tEXdIGcEooSFjYeiEvo5ORKTXxdo0ZGanmNlLZrbSzH6Q53kzs38Pn68zsyNiC+aRH8P0IXDdgbDtHYwwCcCOUUJaRkJEEii2RGBmaeCXwKnAwcC5ZnZwzmWnAgeEP5cAv4orHp7obCioaZSQiCRanE1DRwMr3f01ADObB5wJrMi65kzgt+7uwF/NbC8z28fd3+6xKKYP6fr5Q6fCWfHlHxGRYhdn09AY4K2s4/rwXHevwcwuMbMlZrZkw4YN3Yvi2O8AtM8V6GDQGGhs6N77iYiUmDgTgeU5l1seR7kGd5/t7jXuXjNixIjuRfHZf+3kg1Iw5nAtKCciiRdnIqgH9ss6Hgus2YVr4qMkICISax/BYuAAMxsPrAamAuflXHMfcHnYf3AM8H6P9g9kTH+/x99SRKRUxJYI3L3FzC4HHgbSwFx3f97MLg2fvxl4EDgNWAlsAy6MKx4REckv1gll7v4gQWGffe7mrMcOXBZnDCIi0rXkrjUkIiKAEoGISOIpEYiIJJwSgYhIwpl73jm3RcvMNgBv7OLLhwPv9GA4/YHuORl0z8mwO/f8IXfPOyO33yWC3WFmS9y9pq/j6E2652TQPSdDXPespiERkYRTIhARSbikJYLZfR1AH9A9J4PuORliuedE9RGIiMjOklYjEBGRHEoEIiIJl5hEYGanmNlLZrbSzH7Q1/HEwcz2M7MFZvaCmT1vZrXh+aFm9qiZvRL+uXdfx9qTzCxtZsvM7IHwuNTvdy8zu8vMXgz/rT+WgHv+dvh/ermZ3WFmlaV2z2Y218zWm9nyrHOd3qOZXRmWZy+Z2cm789mJSARmlgZ+CZwKHAyca2YH921UsWgBrnD3g4CPApeF9/kD4DF3PwB4LDwuJbXAC1nHpX6/NwAPufsk4FCCey/ZezazMcC3gBp3/wjBsvZTKb17vg04Jedc3nsMf6+nApPD19wUlnO7JBGJADgaWOnur7l7EzAPOLOPY+px7v62uz8TPm4gKCDGENzrb8LLfgP8XZ8EGAMzGwt8DpiTdbqU73cw8CngFgB3b3L3TZTwPYfKgCozKwMGEuxkWFL37O5/Bt7NOd3ZPZ4JzHP3RndfRbCny9G7+tlJSQRjgLeyjuvDcyXLzKqBw4GngFGZnd/CP0f2YWg97RfA/wHass6V8v1OADYAt4bNYXPMbA9K+J7dfTVwLfAm8DbBToaPUML3nKWze+zRMi0piWDnveuhZMfNmtmewN3ANHff3NfxxMXMTgfWu/vSvo6lF5UBRwC/cvfDga30/yaRLoXt4mcC44F9gT3M7Ct9G1Wf69EyLSmJoB7YL+t4LEHVsuSYWTlBEvi9u98Tnl5nZvuEz+8DrO+r+HrYx4HPm9nrBM19x5vZ7yjd+4Xg/3K9uz8VHt9FkBhK+Z5PBFa5+wZ3bwbuAY6ltO85o7N77NEyLSmJYDFwgJmNN7MBBJ0s9/VxTD3OzIyg7fgFd78+66n7gPPDx+cDf+jt2OLg7le6+1h3ryb4N33c3b9Cid4vgLuvBd4yswPDUycAKyjheyZoEvqomQ0M/4+fQND/Vcr3nNHZPd4HTDWzCjMbDxwAPL3Ln+LuifgBTgNeBl4F/qmv44npHj9BUD2sA54Nf04DhhGMOHgl/HNoX8caw70fBzwQPi7p+wUOA5aE/873Ansn4J5/DLwILAduBypK7Z6BOwj6QJoJvvFf1NU9Av8UlmcvAafuzmdriQkRkYRLStOQiIh0QolARCThlAhERBJOiUBEJOGUCEREEk6JQEQk4ZQIJLHMbJiZPRv+rDWz1eHjLWZ2Uwyfd5uZrTKzS7u45pNmtiJ7KWKRuGkegQhgZtOBLe5+bYyfcRvBpLe7ClxXHV73kbhiEcmmGoFIDjM7LmuTm+lm9hsze8TMXjezs83sGjN7zsweCtd2wsyONLNFZrbUzB7OrA9T4HO+GG608jcz+3Pc9yXSGSUCkcImEux5cCbwO2CBu08BPgA+FyaDG4Fz3P1IYC5wVYT3/RfgZHc/FPh8LJGLRFDW1wGI9AN/dPdmM3uOYHesh8LzzwHVwIHAR4BHgzXRSBOsGVPI/wNuM7P/JlhRU6RPKBGIFNYI4O5tZtbsOzrW2gh+hwx43t0/1p03dfdLzewYgtrGs2Z2mLtv7MnARaJQ05DI7nsJGGFmH4NgTwgzm1zoRWY20d2fcvd/Ad6h4/ryIr1GNQKR3eTuTWZ2DvDvZjaE4PfqF8DzBV46y8wOIKhRPAb8LdZARTqh4aMivUTDR6VYqWlIpPe8D/y00IQy4H6CpiKRXqEagYhIwqlGICKScEoEIiIJp0QgIpJwSgQiIgn3/wGc2lLRc8W4twAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.plot(campsol[:,0],campsol[:,5],label='CAMP')\n",
    "plt.plot(solTchem[:,sp,t_indx],solTchem[:,sp,C_indx],\"*\",label='TChem')\n",
    "plt.xlabel('Time [s]')\n",
    "plt.ylabel('Concentration of B [mol/m3]')\n",
    "plt.legend(loc='best')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
